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Abstract 

The  primary  objective  of  this  study  was  to  develop  an  analytical  methodol¬ 
ogy  based  on  the  Mezm  V'aluc  Analysis  algorithm  that  approximates  the  perfor¬ 
mance  characteristics  of  a  queueing  network  model  (QNM)  containing  a  fork-join 
queue  with  probabilistic  branching.  These  performance  characteristics  are  response 
time,  throughput  and  queue  length  at  each  station.  The  QNM  solved  contains  the 
essential  features  of  the  aircraft  sortie  generation  process.  The  sensitivity  of  the 
method’s  accuracy  to  increases  in  server  utilization  was  determined.  The  compar¬ 
isons  of  the  results  of  the  MVA  heuristic  to  the  outputs  of  the  Logistics  Composite 
Model  (LCOM)  simulation  indicate  that  the  heuristic’s  accuracy  decreased  as  server 
utilization  increases.  When  server  utilization  wa  kept  in  realistic  ranges,  the  results 
of  the  heuristic  for  a  single  fork-join  queue  were  very  accurate.  For  non-maintenance 
stations,  results  were  within  1  to  2  percent  of  the  LCOM  simulation  output.  For 
stations  on  the  fork-join  queue  paths,  heuristic  results  were  within  5  percent  of  the 
LCOM  simulation’s  output  for  that  portion  of  the  network. 


A  MEAN  VALUE  ANALYSIS  HEURISTIC 
FOR  ANALYSIS  OF 
AIRCRAFT  SORTIE  GENERATION 


1.  Introduction 


1  l  Motivation 

Aircraft  sortie  generation  is  at  the  heart  of  the  United  States  Air  Force's  mis¬ 
sion.  The  process  of  aircraft  sortie  generation  includes  all  those  activities  necessary 
for  aircraft  to  fly  the  maximum  number  of  sorties  in  a  specified  period  of  time.  These 
activities  include,  but  are  not  limited  to,  taxiing  and  flying  the  aircraft,  weapons 
loading  or  unloading,  and  repairing  system  malfunctions.  To  successfully  accomplish 
the  mission,  sufllcient  resources,  such  as  equipment  and  personnel,  must  be  available 
to  accomplish  each  activity  in  a  timely  fashion.  On  the  other  hand,  fiscal  contraints 
limit  resource  availability.  The  accurate  determination  of  the  resource  requirements 
has  historically  been  accomplished  through  discrete  event  simulation  of  the  aircraft 
sortie  generation  process.  Analytical  methods  have  not  been  fully  explored.  Perhaps 
one  reason  for  this  is  the  difiiculty  in  analytically  modeling  concurrent  maintenance 
activities. 

1.2  Background 

The  process  of  aircraft  sortie  generation  can  be  thought  of  as  the  operation 
of  a  system,  a  collection  of  entities  (i.c..  aircraft,  peojde,  equipment  and  parts^ 
which  act  and  interact  together.  The  system  description  tan  be  discrctt,  in  that 
the  state  variables  (the  collection  of  variables  necessary  to  describe  the  system  at  a 
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particular  tiriu')  change  only  at  a  countable  (or  finite)  number  of  points  in  time  (20:2). 
1  he  aircraft  sortie  generation  process  can  be  completely  described  by  detailing  the 
number  of  aircraft  in  various  conditions,  such  as  flying,  taxiing,  or  being  repaired. 
I  lie  state  of  the  system  changes  only  when  an  aircraft  moves  from  one  condition  to 
another. 

1.2.1  Modeling  Techniquea.  Because  aircraft  sortie  generation  is  a  discrete 
system  that,  as  will  be  discussed  below,  can  be  accurately  modeled  as  a  queueing  net¬ 
work,  the  most  applicable  tools  for  modeling  and  analyzing  the  system  are  discrete- 
event  simulation  and  queueing  theory  (16:427)  (20:4).  Several  discrete-event  simula¬ 
tions  of  aircraft  sortie  generation  exist.  The  Logistics  Composite  Model  (LOOM)  is 
a  highly  detailed  discrete-event  simulation  that  has  been  validated  and  “designated” 
for  use  by  the  Air  Force  (10).  For  this  reason,  LCOM  is  considered  to  be  the  stan¬ 
dard  by  which  other  models  of  aircraft  sortie  generation  are  judged.  LCOM  and 
other  simulations  of  aircraft  sortie  generation  are  detailed  in  Chapter  2. 

ILi fortunately,  analysis  using  a  simulation  model  can  be  difficult.  The  ex¬ 
haustive  examination  of  alternatives  using  simulation  can  be  very  time  consuming. 
This  situation  is  aggravated  by  the  fact  a  simulation  output  is  random  in  nature 
and  therefore  requires  that  many  iterations  be  run  with  identical  inputs  in  order  to 
establish  a  confidence  interval  for  the  results  (3:87). 

On  the  other  hand,  queueing  theory  can  be  used  to  develop  queueing  network 
models  (QNMs)  that  do  not  require  a  large  amount  of  input  data  and  can  execute 
quickly  on  a  computer  (21:80).  These  models  may  be  analytically  solved  to  obtain 
such  measures  of  performance  as  average  waiting  line  length,  average  waiting  time, 
server  utilization,  and  .system  throughput.  This  capability  allows  the  analyst  to 
quickly  reveal  relationships  between  key  parameters,  evaluate  tradeoffs  and  offer 
fundamental  insights  into  system  performance  (11:16).  For  these  reasons,  QNMs 
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have  been  used  to  analyze  several  aircraft  maintenance-related  problems.  Some 
examples  of  these  applications  are  described  in  Chapter  2. 

The  drawbacks  of  using  QNMs  are  centered  around  the  simplifying  assump 
tions  commonly  made  in  order  to  make  the  model  analytically,  i.e.,  mathematically, 
tractable.  Simplifying  assumptions  are  made  concerning  the  customer  arrival  and 
service  time  distributions,  queueing  discipline,  scheduling  policies,  and  steady-state 
conditions  (8:25).  The  structure  of  the  model  itself  may  prevent  the  use  of  analytical 
methods  for  arriving  at  exact  results,  but  useful  approximate  solutions  can  often  be 
obtained. 

1.3  Problem 

.Aircraft  sortie  generation  is  fairly  independent  of  the  aircraft  type  with  differ¬ 
ences  confined  largely  to  the  maintenance  aspects  of  the  process.  Thus.,  while  the 
number  of  aircraft  systems  and  their  probability  of  failure  is  unique  for  each  aircraft 
type,  the  overall  process  of  sortie  generation  remains  basically  invariant. 

The  process  begins  with  the  aircraft  taxiing  to  the  runway  for  takeoff.  The 
pilot  conducts  system  checks  during  ground  operation  and,  with  a  small  probability, 
some  system  or  systems  may  fail  and  require  maintenance,  thereby  aborting  the 
sortie.  If  the  aircraft  flies,  it  will  do  so  for  a  predetermined  length  of  time.  With 
a  higher  probability  than  that  of  at  the  end-of-runway  check,  some  aircraft  system 
or  systems  may  have  failed  in  flight  and  require  maintenance  upon  landing.  If  no 
maintenance  is  required  upon  landing,  the  aircraft  is  “turned”  (i.e.,  prepared  for 
its  next  sortie)  rearmed  with  weapons,  and  again  taxis  to  the  runway  for  its  next 
sortie.  If  maintenance  is  required  due  to  system  failures  at  the  end  of  the  runway, 
the  aircraft's  weapons  must  be  dearmed  or  removed  before  the  maintenance  can 
begin.  Maintenance,  regardless  of  when  the  failure  occurred,  is  performed  on  all 
malfunctioning  systems  simultaneously.  Each  system  has  a  unique  probability  of 


failure.  Upon  completion  of  all  maintenance,  the  aircraft  is  turned,  rearmed,  and 
taxis  for  its  next  sortie. 

The  entire  process  is  illustrated  in  figure  1.1.  Each  activity  of  the  sortie  gen¬ 
eration  process  can  be  modeled  as  a  station  in  a  queueing  network.  Each  station 
is  numbered,  with  representing  the  probability  that  the  aircraft  transitions  from 
queue  i  to  j.  The  service  times,  s,  ,  and  number  of  servers,  n,  ,  for  each  queue  are  listed 
in  Chapter  4  where  the  network  is  solved.  These  network  parameters  arc  notional; 
only  the  QNM  structure  reflects  an  actual  aircraft  sortie  generation  process. 

Five  aspects  of  the  aircraft  .sortie  generation  process  arc  difficult  to  model  and 
solve  using  queueing  theory.  These  are; 

1.  Concurrent  Maintenance:  Upon  experiencing  a  ground  abort  or  a  sortie  com¬ 
pletion,  an  aircraft  may  have  more  than  one  malfunctioning  system.  These 
systems  could  all  be  repaired  simultaneou.sly  by  personnel  with  separate  main¬ 
tenance  specialties.  This  aspect  of  the  system  is  modeled  with  a  fork-join 
queue.  The  fork-join  queue,  illustrated  in  figure  1.2,  has  multiple  paths  ema¬ 
nating  from  a  fork  node  with  each  path  representing  a  separate  maintenance 
speciality.  The  maintenance  activities  performed  by  the  speciality  are  modeled 
as  a  multiple-server  queue.  The  number  of  servers  represents  the  number  of 
personnel  available  in  that  speciality  and  the  service  time  is  the  average  com¬ 
pletion  time  of  the  maintenance  activities  performed  by  the  speciality.  The 
probabilities,  ,  are  the  conditional  probabilities  that  maintenance  special¬ 
ity  j  must  perform  a  maintenance  activity  given  that  the  aircraft  has  some 
inalfunctioning  system. 

2.  Maintenance  Crew  Sizes:  Personnel  from  a  given  maintenance  specially  may 
be  responsible  for  more  than  one  aircraft  .system.  Repair  of  these  systems 
may  require  differen!  service  times  and  crew  sizes.  This  aspect  of  maintenance 
is  modeled  within  the  fork-join  (jueue  by  averaging  crew  sizes  and  < oinidiu.ion 
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Figure  1.2,  Fork- Join  Queue  for  Maintenance  Activities 
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times  over  all  the  activities  that  each  maintenance  speciality  must  perform.  For 
example,  if  a  maintenance  speciality  has  20  people  available  and  the  average 
crew  size  is  2,  then  the  multiple-server  queue  modeling  that  speciality  htis  10 
servers  and  each  server  is  considered  two  people.  An  alternative  method  of 
addressing  this  issue  is  known  as  the  machine-operator  interference  problem 
and  is  discussed  next. 

3.  Machine-Operator  Interference:  Some  maintenance  tasks,  such  as  de-arming 
aircraft  weapons  before  ground-abort  maintenance  and  rearming  an  aircraft 
prior  to  flight,  occur  at  different  points  in  the  network  but  utilize  the  same 
maintenance  manpower  specialty.  This  linkage  of  resources  is  referred  to  in 
the  literature  of  queueing  network  theory  as  machine-operator  interference. 
In  this  study,  the  issue  is  handled  by  treating  the  different  tasks  as  separate 
queues  with  their  own  resource  pool.  However,  machine-operator  interference 
in  automatic  assembly  systems  with  some  similarities  to  the  aircraft  sortie 
generation  process  h2i8  been  sucessfully  modeled  (7:275)  (17:93). 

4.  Batch  Arrivals:  Some  combat  aircraft  rarely  fly  alone,  but  rather  in  groups 
(flights)  of  two  or  four.  The  simultaneous  arrival  of  multiple  customers  to  a 
fork-join  queue  is  an  aspect  of  the  theory  that  haa  yet  to  be  addressed  in  the 
literature  and  is,  therefore,  not  modeled  in  this  study. 

5.  Mission  Scheduling:  Aircraft  normally  do  not  fly  immediately  after  mainte¬ 
nance  is  complete,  but  instead  fly  according  to  a  piredeterrnined  schedule.  The 
time  aircraft  spend  waiting  for  a  scheduled  take-off  time  reduces  the  flow  of 
aircraft  through  the  network.  Having  aircraft  fly  sorties  as  soon  as  they  can 
creates  a  more  demanding  sortie  generation  scenario.  Maintenance  personnel 
that  can  support  a  “fly  w'hen  ready”  scenario  can  easily  support  a  scenario 
with  delays  due  to  scheduling.  Thus,  mission  scheduling  is  not  modeled  in  this 
study. 
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1.4  Objective 

The  primary  objective  of  this  study  is  to  develop  an  analytical  methodology 
based  on  the  Mean  Value  Analysis  algorithm  that  approximates  the  performance 
characteristics  of  a  queueing  network  model  containing  a  fork-join  queue  with  prob¬ 
abilistic  branching.  These  performance  characteristics  are  response  time,  through¬ 
put  and  queue  length  at  each  queue.  The  QNM  to  be  solved  contains  the  essential 
features  of  the  aircraft  sortie  generation  process.  The  sensitivity  of  the  method’s 
accuracy  to  increases  in  server  utilization  will  be  determined. 

1.4.1  Approach.  The  objectives  of  the  study  are  accomplished  through  the 
following  steps. 

1 .  Develop  a  solution  method  based  on  the  Mean  Value  Analysis  algorithm  that 
approximates  the  solution  of  an  analytical  queueing  network  model  containing 
a  fork-join  queue  with  probabilistic  branching. 

2.  Solve  the  QNM  in  figure  1.1  with  the  fork-node  queue  in  figure  1.2  representing 
maintenance. 

3.  Develop  an  LCOM-version  of  the  same  network.  This  network  will  be  simulated 
in  LCOM  with  a  sufficient  number  of  iterations  to  guarantee  a  reasonably  tight 
confidence  interval  for  the  simulation’s  outputs. 

4.  Compare  the  analytical  solution  of  the  QNM  to  the  LCOM  simulation  output. 

5.  Repeat  steps  2  thru  4  for  the  QNM  with  specific  network  parameters  changed 
such  that  server  utilizations  are  increased. 

1.4.2  Scope.  Aircraft  maintenance  is  divided  into  two  categories  depending 
upon  where  the  maintenance  is  performed.  These  two  categories  are  on- equipment 
and  off-equipment  maintenance,  where  the  term  equipment  refers  to  the  aircraft. 
On-equipment  maintenance  entails  those  tasks  where  the  malfunctioning  aircraft 
system  is  repaired  without  leaving  the  aircraft.  Examples  include  avionics  systems 
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troubleshooting,  component  removal  and  replacement  and  airframe  repair.  Off- 
equipment  maintenance  is  performed  on  systems  that  have  been  removed  from  the 
aircraft  and  taken  to  a  workcenter  for  repair.  The  most  common  instance  of  off- 
equipment  maintenance  is  the  repair  of  an  avionics  system  in  the  corresponding 
workcenter.  This  study  will  address  on-equipment  aircraft  maintenance  only. 
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II.  Previous  Work 


2.1  Introduction 

Previous  work  in  the  modeling  of  aircraft  sortie  generation,  regardless  of  the 
purpose  of  the  modeling,  has  primarily  employed  simulation  methods.  However, 
some  examples  of  analytical  queueing  models  also  exist. 

2.2  Simulation  Models 

Most  military  studies  of  aircraft  sortie  generation  have  relied  on  simulation  as 
the  principal  evaluation  tool  (24).  This  is  primarily  due  to  the  complexity  of  the 
maintenance  systems  that  comprise  a  large  part  of  the  process  (11:2).  For  example, 
the  repair  of  a  single  malfunctioning  aircraft  system  may  require  the  combined  efforts 
of  several  different  maintenance  specialities  working  in  sequence  or  even  simultane¬ 
ously.  Each  of  these  maintenance  specialities  has  a  different  size  manpower  pool 
from  which  to  draw,  a  unique  set  of  aircraft  systems  whose  repair  are  its  responsibil¬ 
ity,  and  its  own  constantly  changing  backlog  of  aircreift  malfunctions  waiting  to  be 
repaired.  Simulation  offers  the  unique  ability  to  capture  the  complexities  of  systems 
that  cannot  be  otherwise  analyzed.  Some  recent  simulations  of  the  aircraft  sortie 
generation  process  are  discussed  below. 

2.2.1  The  Logistics  Composite  Model.  The  Logistics  Composite  Model  (LCOM) 
is  the  aircraft  sortie  generation  model  that  has  been  approved  for  use  in  the  Air 
Force.  It  is  a  discrete  event  simulation  written  in  Simscript  11.5.  LCOM's  original 
purpose,  when  addressed  by  the  Rand  Corporation  in  the  late  1960s,  was  to  provide 
an  analysis  tool  to  relate  base-level  logistics  re.sources  (the  model’s  inputs)  with  each 
other  and  with  sortie  generation  rates  (the  model’s  output)  (14:2).  The  Air  Force 
now  uses  the  model  primarily  to  determine  maintenance  manpower  requirements 
(12:1-1)  (27:1)  (5:3). 
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Maintenance  resource  levels,  i.e.,  spare  parts,  people  f  i.  cl  equipment,  are  estab¬ 
lished  by  the  analyst.  To  obtain  the  minimum  mix  of  resources  required  to  sustain 
the  desired  sortie  rate,  resources  are  constrained  until  the  sortie  rate  is  affected 
(5:10).  However,  because  LCOM’s  results  inherently  contain  the  randomness  of  the 
Monte  Carlo  simulation  process,  the  optimal  solution  is  inexact. 

LCOM’s  size  and  complexity  limit  its  usefulness.  Data  input  is  difficult  and 
time  consuming;  for  a  detailed  weapon  system  study,  an  LCOM  data  base  can  run  to 
several  thousand  lines  of  code  (14:75)  (27:2).  LCOM  runs  are  also  time  consuming; 
a  single  run  from  a  large  database  can  take  over  three  hours.  Since  LCOM  is  run 
iteratively,  with  the  analyst  varying  inputs  in  order  to  determine  changes  in  output, 
LCOM  studies  can  be  several  months  in  duration. 

LCOM  continues  to  be  updated  and  refined  in  response  to  the  demands  of 
long  term  use.  The  model  has  been  made  increasingly  user-  friendly  and,  if  run  on 
modern  computers,  is  relatively  fast  running  (5:14).  LCOM  is  used  by  the  United 
States  Air  Force’s  Air  Combat  Command,  Air  Mobility  Command  and  Air  Material 
Command  (29). 

2.2.S  LMI  Sortie  Generation  Model.  In  the  late  1970s,  The  Logistics  Man¬ 
agement  Institute  (LMI),  working  for  the  Assistant  Secretary  of  Defense  for  Man¬ 
power,  Reserve  Affairs,  and  Logistics,  developed  the  Sortie  Generation  Model  (SGM). 
The  model  Wtis  developed  to  estimate  the  effect  of  varying  levels  of  certain  classes  of 
resources  on  the  readiness  of  tactical  fighter  forces.  The  SGM  system  consists  of  the 
model  itself  and  a  complex  system  of  supporting  software  that  produces  informa¬ 
tion  files  necessary  to  run  the  SGM  such  as  spares  availability  and  the  maintenance 
environment  (1:6-7). 

Scenario  characteristics  are  specified  interactively  by  the  user.  Characteristics 
such  eis  the  number  of  waves  or  flying  periods  per  day,  the  attrition  rate,  the  number 
of  aircraft  in  reserve,  and  the  ground  abort  rate  can  be  set  to  vary  daily.  Other 
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characteristics  such  as  the  sortie  length,  the  probability  that  a  returning  aircraft 
requires  maintenance  and  the  first  and  last  takeoff  times  of  the  day  remain  constant 
throughout  the  scenario.  The  SGM,  a  hybrid  analytic/simulation,  then  estimates 
maximal  sortie-generation  capability  across  time  as  a  function  of  its  input  files  and 
scenario  characteristics  (1:8). 

SGM  accomplishes  this  estimation  by  tracking  the  number  of  aircraft  in  each 
of  five  states:  1)  Mission-capable;  2)  Maintenance;  3)  Not  mission-capable,  supply; 
4)  Combat  loss;  5)  Reserve.  There  are  eight  processes  by  which  aircraft  move  from 
state  to  state.  These  are  determined  by  the  input  files  and  scenario  characteristics 
and  include  ground  aborts,  breakage,  repair,  attrition,  and  parts  demands.  Aircraft 
availability  is  assessed  against  the  daily  flying  schedule  (2:1-3). 

Aircraft  system  malfunctions  in  the  SGM  are  determined  by  a  series  of  random 
draws  from  binomial  distributions.  An  initial  draw  is  made  to  first  determine  if  the 
aircraft  has  any  breakage.  .Additional  draws  are  made  to  determine  what  mainte¬ 
nance  workcenters  have  malfunctioning  equipment  that  require  repair  (2:2-3).  The 
malfunctioning  equipment  is  not  specified.  Maintenance  is  modeled  as  a  multiple- 
server  queuing  process  with  aircraft  tracking  accomplished  via  attribute  marking 
(2:2-7). 

S.S.3  TSAR.  TSAR  is  a  complex  Monte  Carlo  simulation,  written  in  FOR¬ 
TRAN,  of  a  system  of  interdependent  theater  airbases.  TSAR  was  developed  by 
the  Rand  Corporation  in  the  early  1980s  for  analyzing  the  interrelations  among  on- 
bcise  resources  and  the  capability  of  the  airbases  to  generate  aircraft  sorties  in  a 
wartime  environment.  On-equipment  maintenance  tasks,  parts  and  equipment  re¬ 
pair  jobs,  munitions  assembly  and  facility  repair  tasks  are  simulated  for  each  of 
several  airbases.  Asset  accounting  for  each  of  11  classes  of  resources,  with  hundreds 
of  resources  in  each  class,  permits  assessment  of  a  broad  range  of  policy  options  that 
could  improve  the  efficiency  of  resource  utilization  on  a  theater-wide  basis  (13:v). 


TSAR  can  also  model  the  effects  of  conventional  or  chemical  airbase  attacks,  the 
effects  of  having  personnel  operate  in  individual  chemical  protection  equipment,  and 
the  efforts  to  reconstitute  the  airbase  following  an  attack  (23:4). 

TSAR  is  a  large  and  complex  model  requiring  a  lot  of  upkeep.  There  are  no 
data  preprocessors  for  TSAR  that  are  able  to  take  pre-collected  data  and  transform 
it  into  TSAR  format;  thus,  the  process  of  data  preparation  and  input  can  be  very 
time  consuming  (23:6).  Just  as  in  the  case  of  LCOM,  TSAR’s  results  inherently 
contain  the  randomness  of  the  Monte  Carlo  simulation  process;  several  iterations  of 
the  model  must  be  run  in  order  to  establish  a  confidence  interval  for  the  output  data. 

2.2.4  Dyna-Sim.  Dyna-Sim  is  a  discrete-event  simulation  developed  by  the 
Rand  Corporation  in  the  mid-1980s.  The  model,  written  in  Simscript  II. 5,  was  built 
to  simulate  the  queueing  of  aircraft  maintenance  jobs  for  automatic  test  eijuipment 
(ATE)  in  order  to  quantify  the  importance  of  ATE  in  the  repair  cycle  (22:iii).  The 
simulated  repair  system  is  a  multi-server,  multi-job-class  queue  with  time-varying 
arrival  rates.  The  simulated  queue  is  multi-server  because  more  than  one  ATE 
is  normally  available  for  repair  actions.  The  multi-job-class  «ispect  of  the  model 
captures  the  fact  that  the  ATE  is  used  to  repair  more  than  one  type  of  malfunctioning 
system.  Jobs  are  selected  from  the  maintenance  queue  using  four  different  priority 
rules  (22:6). 

While  the  scope  of  the  Dyna-Sim  model  is  limited  to  only  the  .ATE  issue, 
it  could  be  easily  adapted  to  model  any  maintenance  repair  activity  with  multiple 
servers  and  multiple  jcb-classes.  However,  the  analyst  must  specify  the  various  job 
arrival  rates,  i.e.,  the  system  does  not  model  the  closed  network  aspects  of  aircraft 
sortie  generation.  The  analyst  specifies  when  a.  ival  rales  change  rather  than  have 
these  changes  dependent  upon  repair  frequencies,  service  times,  and  queue  lengths 
(22:3). 
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8.5  Analytical  Models 


Analytical  models  have  been  little  used  in  the  modeling  of  aircraft  sortie  gen¬ 
eration.  Their  use  has  been  primarily  confined  to  small  aspects  of  the  larger  system 
such  as  specific  manpower  or  equipment  issues.  Some  examples  of  analytical  methods 
applied  to  the  aircraft  sortie  generation  process  are  discussed  below. 

2.5.1  Fleet  Maintenance.  Boling  and  Hillier  used  queueing  theory  to  analyze 
fleet  maintenance  systems,  where  the  fleet  consists  of  trucks,  aircraft  or  cars  on  a 
scheduled  maintenance  program  as  well  as  emergency  repair  basis  (4:1).  They  devel¬ 
oped  a  general  queueing  network  model  for  the  design  of  fleet  maintenance  systems 
involving  several  crews  working  in  sequence.  Emphasis  was  placed  on  optimizing  the 
system  for  scheduled  periodic  maintenance  activities.  The  optimal  system  configu¬ 
ration  was  determined  as  a  function  of  the  number  of  crews,  cost  of  crew  activities 
and  cost  of  a  fleet  unit’s  idle  time  (4:2).  In  addition  to  the  assumption  of  sequential 
tasks,  the  method  presented  is  also  based  on  the  assumption  that  the  expected  ser¬ 
vice  times  for  the  sequential  tasks  are  equal.  These  assumptions  could  generally  not 
be  made  in  developing  a  model  of  a  sortie  generation  process  for  military  aircraft. 

2.5.2  B-52H/KC~155  Maintenance  Model.  A  1989  study,  conducted  by  the 
.Air  Force  Logistics  Command,  examined  the  maintenance  manpower,  spares  cost 
and  repair  cost  impacts  on  going  from  three  to  two  levels  of  maintenance  on  avionics 
spares  for  the  B-52H  Bomber  and  the  KC-135  Tanker.  Queueing  theory  was  used 
to  show  that  consolidating  several  bcise  maintenance  shops  into  a  single  repair  fa¬ 
cility  would  reduce  the  maintenance  manpower  needed  to  give  the  same  or  better 
service  (25:1).  This  was  achieved  by  bounding  the  number  of  servers  required  in  the 
repair  facility  from  above  ari'  ^elow.  The  upper  bound  was  determined  by  solving 
an  M/M/s  queue  for  the  iiunimum  number  of  servers  required  to  ensure  that  the 
expected  waiting  time  under  consolidation  was  less  than  or  equal  to  the  base  ex¬ 
pected  waiting  time.  The  lower  bound  was  determined  by  combining  servers  into 


a  consolidated  servei  with  a  service  rate  equal  to  the  base  service  rate  multiplied 
by  the  number  of  servers  combined.  The  expected  waiting  time  was  held  constant 
and  the  resulting  M/M/1  queue  was  solved  for  the  fraction  of  servers  that  must  be 
retained  under  consolidation  (25:45-49).  While  a  heuristic  argument  for  the  accep¬ 
tance  of  these  bounds  was  presented,  no  formal  proof  or  validation  of  the  bounds 
was  provided. 

2.3.3  Optimal  Specialization.  Dietz  and  Rosenshine  developed  an  analytical 
method  for  determining  an  optimal  specialization  strategy  for  a  maintenance  man¬ 
power  force.  The  model  assumed  that  maintenance  tasks  were  generated  by  a  system 
of  identical  machines,  such  as  aircraft,  which  experience  reindom  malfunctions  and 
require  periodic  service.  As  a  part  of  the  overall  method,  the  impact  of  alternative 
manpower  structures  on  system  performance  was  evaluated  using  a  queueing  net¬ 
work  model  (11:2).  The  method  requires  the  simultaneous  solution  of  the  network’s 
global  balance  equations,  thus  limiting  the  applicability  of  the  approach  to  simple 
systems  or  highly  aggregated  analysis  of  complex  systems  (11:16). 

2.3.4  SUMMA.  The  Small  Unit  Maintenance  Manpower  Analyses  (SUMMA) 
method  wcis  developed  in  the  late  198Us  for  the  Air  Force  Human  Resource  Lab.  The 
method  was  originally  intended  to  be  used  for  gauging  the  impacts  of  airbase  dis¬ 
persal  in  wartime  on  manpower  requirements  and  unit  performance  under  proposed 
alternative  arrangements  of  maintenance  job/task  skills  (6:i).  The  SUMMA  method 
grew  into  a  comprehensive  system  that  would  consider  maintenaTice,  manpower,  per¬ 
sonnel  and  training  (MPT),  and  cost  factors  to  help  the  Air  Force  improve  ii  main¬ 
tenance  task  allocation/spccialty  structure.  The  system  uses  maintenance  task  and 
MPT  data  for  specific  operational  and  maintenance  scenarios  in  a  decision  support 
system  (DSS)  to  expose  the  tradeoffs  and  implications  of  occupational  specialty  con- 
soliddtion  and/or  other  maintenance  task/specialty  changes  (30:6).  The  SUMM.\ 
DSS  analyzes  the  maintenance  process  in  terms  of  task  completion  times.  Using 


um‘  raclhod  of  Lagrangian  Multipliers,  the  number  of  personnel  in  each  speciality  is 
n’iinimizcd  subject  to  maximum  aircraft  downtime  constraints.  Thus  a  task/specialty 
a.ss’.gnmcnt  that  achieves  the  desired  sortie  rate  using  the  minimum  manpower  while 
meeting  these  user-specified  constraints  is  derived  (30:99-110).  The  operational  ca¬ 
pability  of  the  task  allocation  is  then  verified  using  LCOM. 

The  SUMMA  method  contains  a  major  assumption  that  does  not  reflect  real- 
world  practice.  In  the  SUMMA  task  allocation  model  it  is  assumed  that  all  tasks  on 
an  aircraft  are  performed  sequentially  (30:104).  This  assumption  is  enforced  even 
if  the  aircraft  requires  maintenance  on  separate  systems  maintained  by  separate 
specialities  that  could  work  on  the  aircraft  simultaneously.  In  actual  practice,  these 
systems  would  be  repaired  at  the  same  time.  The  SUMMA  model  itself  relies  on 
LCOM  to  evaluate  and  verify  the  task  allocation  derived  by  the  SUMMA  method. 
The  SUMMA  method  produces  “...ttisk/AFS  allocations. ..judged  to  have  flaws  that 
cannot  be  accepted”,  thus  the  method  “...indeed  expects  the  user  to  revise  and  refine 
hi:-i  solution  a  number  of  times”  (30:21). 


III.  Methodology 


3.1  Introduction 

The  method  used  in  this  study  to  analyze  the  aircraft  sortie  generation  process 
is  a  heuristic  based  upon  the  Mean  Value  Analysis  (MVA)  algorithm.  The  MVA 
algorithm,  which  yields  exact  results  for  closed  networks  that  have  a  product  form 
solution,  can  be  modified  to  produce  approximate  results  for  networks  which  do  not 
have  product  form  solutions,  such  as  a  sortie  generation  model  (18:33).  A  discussion 
of  tne  accuracy  of  the  heuristic  is  postponed  until  Chapter  4  when  output  of  the 
LCOM  simulation  of  the  network  is  compared  to  the  results  attained  from  the  MVA 
heuristic.  The  terms  ‘‘queue”  and  “station”  are  used  interchangably  throughout  this 
chapter. 

3.2  The  MVA  Algorithm 

The  MVA  algorithm  for  closed  i  etworks  with  product  form  solutions  is  a  recent 
development  in  queueing  theory,  first  published  in  1980  (19).  It  allows  us  to  analyze 
the  steady-state  behavior  of  a  closed  network  that  exhibits  flow  balance  and  one-step 
behavior.  A  network  exhibits  flow  balance  if  the  total  number  of  arrivals  to  a  station 
during  a  period  T  equals  the  total  number  of  departures  from  the  station.  One-step 
behavior  means  that  arrivals  do  not  coincide  with  departures  and  at  any  instant  only 
one  arrival  or  one  departure  can  occur  (18:89). 

A  closed  network  consists  of  A'  customers  distributed  at  K  interconnected 
queues  (stations).  The  K queues  ma  be  single-server  (load-independent)  or  multiple- 
server  (load-dependent).  The  service  rate  at  a  station  i  is  denoted  fi,  while  the  service 
time,  the  inverse  of  the  service  rate,  i  lenoted  s,.  P,j  is  the  probability  that  a  cus¬ 
tomer  completing  service  at  station  i  will  then  proceed  to  station  j.  The  response 
time  of  a  .station  i.  Hi,  measures  the  combined  waiting  and  service  time  for  a  cus¬ 
tomer.  The  queue  length,  Q,,  is  the  number  of  cutomcr.s  waiting  for  service  and 
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being  served  at  station  i.  The  utilization  of  station  i,  Ui,  has  two  interpretations.  If 
station  i  is  a  single-server  queue,  it  is  the  average  percentage  of  time  the  server  is 
busy  or,  if  the  station  is  a  multiple-server  queue,  it  is  the  average  number  of  busy 
servers. 

The  MVA  algorithm  yields  the  mean  values  of  the  network  performance  charac¬ 
teristics  of  response  time,  queue  length,  throughput,  and  utilization  for  each  station 
in  the  network.  Mean  values  are  normally  the  most  desired  outputs  of  a  model 
(15:122).  Simply  put,  the  algorithm  determines  the  mean  response  time  of  a  sta¬ 
tion  i  in  a  network  to  be  the  sum  of  the  mean  service  times  and  the  mean  waiting 
times  a  customer  experiences  before  being  served.  This  relationship  is  based  on  two 
fundamental  results;  the  .Arrival  Theorem  and  the  Marginal  Local  Balance  Theorem. 

The  Arrival  Theorem  tells  us  that  an  arriving  customer  sees  a  random  ob¬ 
server’s  distribution  with  the  arriving  customer  removed  from  the  network.  This  can 
be  stated  as: 

PA,(n|A')=  P,(nlAf- 1)  (3.1) 


where  P/l,(nlAf)  is  the  probability  that,  in  a  closed  network  containing  A  customers, 
an  arriving  customer  to  station  i  finds  n  customers  already  there  and  P,(n|A'  —  1)  is 
the  probability  that  a  random  observer  sees  n  customers  at  station  i  when  there  are 
N  -  1  customers  in  the  network  (18:97). 

From  the  Arrival  Theorem  we  get  the  Marginal  Local  Balance  Theorem,  which 
can  be  interpreted  as  a  detailed  balance  between  two  adjacent  states  for  station  i  as 
shown  in  figure  3.1. 


The  rate  of  departure  from  state  (n|A')  is  /i,(n),  whereas  the  rate  of  entry  into 
it  from  state  (n  — 1|A  —1)  is  the  throughput  A, (A’),  yielding  a  local  balance  equation 
of 
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Figure  3.1.  Marginal  Local  Balance  for  Station  i 


We  also  not^  that 


/i,(n)  =  HiCi(n)  = 


C,{n) 


where  is  the  basic  service  rate  and  s°  is  the  basic  service  time  for  a  server  at 
station  i  and  C,(n)  is  the  number  of  servers  busy  when  n  customers  are  at  station  i 
(18:232).  From  this  we  get  a  statement  of  the  Marginal  Local  Balance  Theorem: 


-l|iV-l)  =  M^P,(„-l|,V_l)  (3.2) 


C,(n) 


where  the  last  equation  makes  use  of  the  Utilization  Law  (18:94,102), 


Ui{n)  =  s- Ai(n) 


This  theorem  can  be  used  to  derive  several  other  relationships  that  are  useful  in 
computing  performance  measures  recursively.  For  example, 

g,(A')  =  f:  nP.(n|A)  =  -  1|.\'  -  1)  (3.3) 
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Recalling  Little’s  Law  for  a  station  i  when  n  customers  are  present, 


Q,(n)  =  Xi{n)Ri{Ti) 


and  again  making  use  of  the  Utilization  Law  (18:232),  we  substitute  Equation  3.3 
into  Little’s  Law  and  solve  for  the  response  time  to  yield 

RAN)  =  s:'£-^P.{n-l\N-l)  (3.4) 

If  station  i  has  only  one  server  then  C,{n)  =  1  for  all  n  and  the  same  substitution 
vields, 

R,{N)  =  s,  Y.  -  1|A^  -  1)  =  a.[l  +  Q.(A'  -  1)]  (3.5) 

ri=l 

where  Q,(0)  =  0.  Note  that  if  station  i  has  more  servers  than  potential  customers, 
the  station  is  referred  to  as  a  delay  station.  In  this  situation,  Equation  3.5  yields 
Ri{n)  =  Si. 

Equation  3.5  is  known  as  the  mean  value  theorem.  Equations  3.4  and  3.5  are 
at  the  heart  of  the  MVA  Algorithm.  By  relating  the  response  time  of  a  station 
when  n  customers  are  present  to  the  length  of  the  queue  at  the  station  when  n  -  1 
are  present,  we  are  able  to  iteratively  determine  each  station’s  steady  state  queue 
length,  response  time,  and  utilization.  In  order  to  make  these  determinations,  it  is 
necessary  to  calculate  the  average  cycle  time  for  a  customer  when  n  customers  are 
in  the  system.  A  customer  at  station  1  is  chosen  as  reference.  The  time  between  two 
departures  by  the  sa  ne  customer  from  station  1  is; 

=  (3.6) 


where  —  is  the  mean  number  of  visits  customers  make  to  station  i  for  every  visit  to 


station  1.  The  values  u,  are  referred  to  as  visit  ratios.  Visit  ratios  are  determined 
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by  solving  the  system  of  equations  v  —  vP  where  P  is  the  matrix  of  transition 
probabilities  for  the  network.  In  order  to  solve  this  system,  one  of  the  i>,  can  be  set 
to  an  arbitrary  value  (normally  Ui  =  1). 

The  throughput  for  station  1,  when  n  customers  are  in  the  network,  is 

A,(n)=j^  (3.7) 

The  ratio  —  determines  the  amount  of  station  1  throughput  at  station  i.  Thus,  Ai(n) 
is  used  ^  determine  the  queue  length  and  utilization  at  each  station  as  follows: 


Q,{n)  =  R,{n)\,{n)^ 


(3.8) 


and 

ih{n)  =  sMn)-  (3.9) 

Vl 

If  station  i  has  a  single  server,  then  the  result  of  Equation  3.8  can  now  be 
applied  to  Equation  3.5  to  begin  the  next  iteration  of  the  MVA  algorithm.  If  station  i 
h2is  multiple  servers,  the  result  of  Equation  3.9  and  the  Marginal  Local  Balance 
Theorem  can  be  used  to  find  P,(n  —  l|jV  —  1)  as  follows: 


P,{k\n) 


U,{n)Pi{k  -  l|n  -  1) 
C,{k) 


(3.10) 


and 

ft(0|r.)  =  1  -  y:  P,(tW  (3.11) 

*  =  1 

These  two  results  can  now  be  used  in  Equation  3.4  to  begin  the  next  iteration  of  the 
MVA  Algorithm. 

In  summary,  the  MV'A  algorilhni  begins  by  noting  that  Q,{0)  =  0  so  that 
i?,(l)  =  Si.  Equations  3.6  thru  3.9  are  calcualted  t<>  determine  the  values  of  the  mean 
performance  characteristics  of  each  station  Equations  3.10  and  3.11  are  calculated 
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for  multiple  server  stations.  The  next  iteration  of  the  MVA  algorithm,  i.e.,  when  the 
number  of  customers  in  the  networks  is  n  =  2,  uses  the  information  of  the  previous 
iteration  as  detailed  above.  Iterations  continue,  updating  the  values  of  the  mean 
performance  characteristics  at  each  station  with  each  iteration,  until  n  =  A^  The 
MVA  algorithm  is  then  complete  and  the  final  values  foi  each  station’s  response 
time,  throughput,  queue  length,  and  utilization  represent  the  mean  performance 
characteristics  of  the  network  under  study. 

3.3  An  MVA  Heuristic  for  Fork- Join  Queues 

In  order  to  model  maintenance  activities  occurring  simultaneously  within  the 
aircraft  sortie  generation  process,  it  is  necessary  to  add  a  fork-join  queue  to  the  closed 
queueing  network.  The  new  network  no  longer  possesses  a  product  form  solution, 
so  the  MVA  algorithm  can  not  be  applied  (19:1048).  Recently,  however,  a  heuristic 
based  on  the  MVA  was  developed  by  Rao  and  Suri  to  deal  with  this  case  (26).  The 
following  discussion  summarizes  their  work. 

The  fork-join  queue  in  the  network  studied  by  Rao  and  Suri  consists  of  an 
independent  single  server  queue  on  each  path  ememating  from  a  fork  node.  Service 
times  of  the  queues  need  not  be  identical.  Each  customer  in  the  network  takes 
each  path  from  the  fork  node.  A  customer  cannot  leave  the  join  node  until  the 
service  activities  of  each  path  are  completed  for  that  customer.  That  is,  a  parent 
customer  breaks  into  statistically  identical  sibling  customers  at  the  fork  node,  one 
sibling  per  path.  Each  sibling  customer  finishing  service  must  wait  until  all  other 
sibling  customers  have  completed  their  service  on  their  fork  path.  Once  the  sibling- 
customers  have  all  completed  service,  they  are  re-joined  into  the  parent  customer  at 
the  join  node.  A  typical  fork-join  queue  is  illustrated  in  figure  3.2  below.  Stations  1 
to  K  represent  the  I\  single-server  queues  on  the  paths  of  the  fork-join  node  while 
station  K  -V  \  is  a  single-server  queue  that  comprises  the  remainder  of  the  network. 
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The  Hi  represent  the  queueing  areas  and  the  are  the  waiting  areas  after  service 
for  the  A'  paths. 


Figure  3.2.  Closed  Netw'ork  with  Fork-Join  Queue 

The  MVA  heuristic  for  fork-join  queues  makes  two  key  assumptions  (26:8). 
These  two  assumptions  are  exact  when  only  one  customer  is  in  the  network.  When 
more  than  one  customer  is  present  these  assumptions  become  approximations,  hence 
the  MVA  algorithm  becomes  a  heuristic.  The  first  assumption  is  a  heuristic  appli¬ 
cation  of  the  Arrival  Theorem.  It  supplies  the  pivotal  equation  used  in  the  iterative 
calculations  of  the  mean  peformance  characteristics  of  the  network. 

(Al)  An  arriving  customer  to  a  path  of  the  fork-join  queue  sees  a  random  observer’s 
distribution  with  itself  removed  from  the  network  (that  is,  removed  from  each 
of  the  fork-join  paths). 

This  can  be  expressed  as: 

F.4.(n)  =  F.(n-l)  (3.12) 
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where  n  is  the  population  vector  for  the  fork-join  queue,  with  n*  =  n  Vfc  and  n  -  1 
is  the  population  vector  for  the  fork-join  queue  with  =  n  —  1  Vfc.  As  with  the 
MVA  algorithm  this  leads  to 


R,{n)  =  s,[l  -f  Q{n  -  1)] 


(3.13) 


where  Qi(0)  =  0.  Because  Rao  and  Suri  considered  only  the  case  where  all  sta¬ 
tions,  station  A'  -f  1  as  well  as  the  fork-join  paths,  are  single-server  queues,  a  load- 
dependent  version  of  Equation  3.13  is  not  needed  (26:1).  Their  heuristic  will  be 
expanded  to  cover  multiple-server  queues  on  the  fork-join  paths  in  the  next  section. 

The  second  assumption  is  critical  to  determing  the  length  of  time  sibling  cus¬ 
tomers  must  wait  until  they  can  be  re-joined  into  the  parent  customer  at  the  join 
node. 


(A2)  The  response  time  distribution  of  a  customer  at  station  k  is  exponential  and 
independent  of  the  response  time  of  its  sibling  customers  on  the  other  paths. 


Using  fork-join  queue  path  1  as  reference,  we  cein  determine  how  long  sibling 
customer  1  has  to  wait  until  all  of  its  sibling  customers  complete  service  and  then  may 
be  rejoined  into  the  parent  customer.  The  probability  that  sibling  customer  1  must 
wait  for  sibling  customer  2  is  merely  the  probability  that  sibling  customer  1  finishes 
service  before  sibling  customer  2.  Using  the  well  known  formula  for  calculating  this 
probability  when  the  service  rates  are  exponential,  we  have 


Pi(n) 


mi(n)  _ 

mi(n)  -!-  m2(n) 


(3.14) 


where  m,  is  the  service  rate  of  fork-join  queue  path  i.  Assumption  (A2)  allows  us 
to  apply  the  memoryless  property  of  the  exponential  distribution  to  determine  that 
the  time  sibling  customer  1  must  wait  for  sibling  customer  2  is  /?2(n). 
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Continuing  to  use  fork-join  queue  path  1  as  reference,  the  probability  that 
sibling  customers  1  and  2  must  wait  for  sibling  customer  3  is 

P3(n)  =  P{h{n)  >  max[<i(n),t2(n)]} 

where  the  ti(n)  are  independent  exponential  random  variables  with  rates  7n,(n).  The 
cumulative  distribution  function  of  V2  =  max[ti(n),  ^2(0)]  is  given  by: 

Fy,{x)  =  (1  -  e-"" <“)')(!  - 


Differentiating  with  respect  to  x,  we  obtain  the  probability  density  function  of  Tj. 
fy,  =  m,(n)e-”“(“)"  +  -  [mi(n)  -f 

Thus, 

P3(n)  =  P{<3(n)  >  >2}  =  /  P{h{n)  >  x}fY^{x)dx 

Jo 

Now,  since  t3(n)  is  exponentially  distributed  with  rate  7713(11),  we  have 

fOO 

P3(n)  =  / 

Jo 

-[mi(n)  +  dx 


so  that 


P3(n) 


mi{n) 

mi(n)  +  m3{n) 


m2(n) 

m2(n)  -I-  m3(n) 


mi(n)  +  m2(n) 
mi(n)  -!-  m2(n)  -f  7773(0) 


(3.15) 


Again  using  Assumption  (A2),  we  determine  the  wailing  time  to  be  /?3(n).  For  each 
of  the  k  paths,  the  probabilities  p,(n)  and  average  waiting  times  /?,(n)  that  sibling 
customers  1  through  i  -  1  must  wait  for  sibling  customer  i,  with  7  =  2,3, .. .  ,k.  are 
calculated  in  a  similar  manner  (26:10-12). 
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We  are  now  able  to  calculate  the  mean  time  a  customer  spends  in  the  system. 
Considering  the  subnetwork  consisting  only  of  fork-join  path  1  and  station  I\  -f  1, 
the  time  between  two  departures  from  station  1  for  the  same  customer  when  n 
customers  are  in  the  subnetwork  is 

K 

Ti{n)  =  Ri{n)  +  +  ^K+i(n)  (3.16) 

1=2 

We  can  now  apply  Little’s  law  to  this  subnetwork  to  get 

A,(n)=5^  (3.17) 

Since  arrivals  to  the  fork-join  queue  paths  are  simultaneous,  the  throughputs  of  all 
the  subnetworks  consisting  of  a  single  fork-join  path  and  station  K  +  1  are  identical. 
That  is, 

Ai(n)  =  A2(n)  =  •  •  •  =  AK(n)  =  A(n)  (3.18) 

where  A(n)  is  the  network  throughput  when  n  customers  are  present.  We  also  note 
from  the  last  two  equations  that 


Ti(n)  =  T2(n)  =  •  •  •  =  Ti,{n)  =  T{n)  (3.19) 

which  justifies  our  use  of  the  subnetwork  consisting  only  of  fork-join  path  1  and 
station  +  1  as  reference. 

Each  station’s  utilization  and  queue  length  can  be  calculated  as 


U,{n)  =  3,A(n) 


(3.20) 


and 


Q,{n)  =  A,(n)A(n) 


(3.21) 
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The  resutlts  of  Equation  3.21  can  now  be  substituted  into  Equation  3.13  to  begin 
the  next  iteration  of  the  MVA  heuristic  for  fork-join  queues. 

Just  as  in  the  MVA  algorithm,  the  MVA  heuristic  for  fork-join  queues  begins 
by  noting  that  Q,(0)  =  0  so  that  /?,(!)  =  s,.  Next,  Ti(l)  is  calculated,  allowing  the 
throughput  when  one  customer  is  in  the  system  to  be  determined.  The  utilization 
and  queue  length  of  each  station  is  then  calculated  using  Equations  3.20  and  3.21. 
This  information  is  used  in  Equation  3.13  in  the  next  iteration  for  n  =  2  customers. 
Iterations  continue  until  n  =  N  when  t)ie  mean  performance  characteristics  of  each 
station  will  have  been  determined. 

3.4  Extension  to  Fork-Join  Queues  with  Probabilistic  Branching 

A  QNM  of  aircraft  sortie  generation  is  much  more  complicated  than  the  basic 
fork-join  QNM  discussed  above.  The  method  used  to  solve  the  aircraft  sortie  gener¬ 
ation  QNM,  while  similar  to  that  of  the  ba.sic  fork-join  QNM,  must  account  for  these 
complications.  For  example,  in  the  above  calculations,  we  assumed  that  a  customer 
took  each  path  of  the  fork-join  queue.  In  an  aircraft  sortie  generation  QNM,  where 
the  fork-join  queue  represents  the  maintenance  activities  that  can  occur  on  the  air¬ 
craft  simultaneously,  branching  at  the  fork  node  is  probabilistic.  That  is,  when  an 
aircrait  has  a  malfunction,  it  may  have  any  combination  of  system  malfunctions. 
Thus,  the  aircraft  “customer”  may  fork  into  sibling  customers  that  take  fork-join 
paths  1,  2,  and  4  or  paths  3  and  4  or  paths  1,  2,  5  and  6  and  so  forth.  We  must 
account  for  this  probabilistic  branching  in  determining  the  queue  lengths  and  server 
utilizations  for  each  fork-join  path. 

The  response  time  of  queueing  stations  outside  of  the  fork-join  queue  are  deter¬ 
mined  using  the  MVA  algorithm  equations  for  single-server  or  multiple-server  ciucues 
as  appropriate.  The  response  time  of  each  fork-join  queue  path  is  calculated  using 
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MVA  algorithm  Equation  3.4  for  load-dependent  stations, 

n=l 

and  not  Equation  3.13  for  load-independent  stations,  as  in  the  network  studied  by 
Raxj  ajid  Suri.  This  is  because  the  basic  fork-join  QNM  and  the  aircraft  sortie 
generation  QNM  differ  in  the  number  of  servers  on  each  path  of  the  fork-join  queue. 
Each  aircraft  maintenance  speciality  will  have  a  pool  of  personnel  to  draw  from. 
Thus,  instead  of  a  single-server  queue  on  each  path,  the  aircraft  sortie  generation 
fork-join  queue  heis  a  multiple-server  queue  on  each  path. 

Due  to  the  probabilistic  branching  of  the  fork-join  queue,  the  cycle  time  of  a 
representative  customer  cannot  be  determined  as  in  the  MVA  heurisitic  developed 
by  Rao  and  Suri.  Recall  that  they  calculated  the  cycle  time  from  a  single  subnetwork 
consisting  of  a  representative  fork-join  queue  path  and  the  rest  of  the  network  outside 
the  fork-join  queue.  We  must  determine  this  cycle  time  for  K  different  subnetworks, 
one  for  each  fork-join  queue  path  and  the  rest  of  the  network  outside  the  fork-join 
queue. 

We  begin  by  deterimining  the  average  response  time  of  the  fork-join  queue 
given  that  a  customer  takes  path  i.  From  the  first  two  terms  of  Equation  3.16  we 
can  determine  the  time  that  a  sibling  customer  spends  on  fork-join  queue  path  i 
Waiting  to  be  served,  being  served,  and  waiting  for  its  other  sibling  customers  to 
finish  their  service.  This  time  is; 

r.(n)  =  /?,(n)  +  (3.22) 

where  j  =  1, 2, . . . ,  A.  As  discussed  in  the  previous  section,  P;(n)  is  the  probability 
that  sibling  customers  1  through  j  -  1  must  wait  for  sibling  customer  j  and  /?j(n), 
the  response  time  on  fork-path  j,  is  how  long  they  must  wait.  The  term  Pp  is  the 
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probability  that  fork-join  path  j  is  taken,  given  that  at  least  one  of  the  paths  are 
taken. 

The  cycle  time  for  station  1  of  a  representative  customer  through  fork-join 
path  i  is 

C7'(n)  =  fl„(n)—  +  Y.  ^k(n)-  (3-23) 

where  /2/j(n)  =  T,(n)  and  —  is  the  visit  ratio  to  the  fork-join  queue.  The  term 
is  the  sum  of  the  response  times  for  the  stations  outside  of  the  tork- 
join  queue,  adjusted  by  their  respective  visit  ratios. 

Before  we  can  determine  the  average  cycle  time  for  the  QNM  we  need  one  more 
cycle  time,  that  for  the  additional  subnetwork  that  bypasses  the  fork-join  queue  and 
consists  only  of  the  stations  outside  of  the  fork-join  queue  that  do  not  involve  the 
activities  of  the  fork-join  queue.  Any  station  i  with  transition  probability  P,  jj  =  1  to 
station  fj.  the  fork-join  queue,  is  not  included  in  this  final  subnetwork  (for  example, 
in  the  network  illustrated  in  figure  1.1,  station  3,  de-arm  weapons,  is  not  included 
in  this  final  cycle  time  calculation  because  including  it  forces  the  inclusion  of  the 
fork-join  queue).  This  final  cycle  time  is 

C7;+,(n)  =  X:fl*(n)-  (3  24) 

k 

the  sum  of  the  response  times  for  the  k  stations  outside  of  the  fork-join  queue  that 
do  not  involve  the  activities  of  the  fork-join  queue,  as  described  above,  adjusted  by 
their  respective  visit  ratios. 

In  order  to  determine  the  QNM’s  cycle  time  we  average  the  cycle  times  calcu¬ 
lated  for  each  subnetwork  weighted  by  their  respect!  e  visit  ratios: 

K  ,Jj  ,,/j 

CTQ^Min)  =  TCT,{tx)P/^—  +  CT,^,{n){l - )  (3.25) 
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Each  CTi  is  adjusted  by  the  fork-join  node’s  visit  ratio, - ,  and  the  probability  that 

fork-join  queue  path  i  is  the  one  path  definitely  taken,  this  being  Similarly, 
CTi+i  is  adjusted  by  the  visit  ratio  for  the  subnetwork  that  bypasses  the  fork-join 


queue,  which  is  (1  -  — To  preclude  the  possibility  of  the  the  term  (1  —  — )  being 


!} 


’^1 


negative,  vi  is  selected  to  be  a  station  in  the  QNM  with  maximal  throughput. 

As  in  Equation  3.7,  we  determine  the  throughput  through  station  1  when  n 
customers  are  in  the  system  to  be 


A,(n) 


n 

CTQNM{n) 


(3.26) 


For  stations  outside  the  fork-join  queue  (non-maintenance  stations),  the  throughput 
is  applied  to  Equations  3.8  and  3.9  to  determine  the  the  mean  queue  lengths  and 
server  utilizations.  However,  the  throughput  may  not  be  directly  applied  to  Equa¬ 
tions  3.20  and  3.21  to  determine  the  utilization  and  queue  length  of  each  fork-join 
path  queue.  We  must  again  account  for  the  probabilistic  branching  in  the  fork-join 
queue.  This  is  accomplished  by  including  terms  for  the  visit  ratio  to  the  fork-join 
queue  and  the  probability  of  each  path  being  taken.  Thus,  for  each  station  i  in  the 
fork-join  queue,  EqiiJ»t,ions  3.20  and  3.21  become: 


U,{n)  =  s,Ai(n)— P/-' 
1^1 


(3.27) 


and 


Q,(n)  =  K.{n)X,{n)—P!’ 

Vi 


(3.28) 


wh  .e  -  is  the  visit  ratio  to  the  fork-join  queue  and 

fork  join  path  i  ij  taken.  In  aircraft  sortie  generation  terms.  —  is  the  number  of 

Vl 


p!‘  is  the  probability  that 


aircraft  that  experience  a  system  malfunction  for  every  aircraft  that  taxis,  while 
is  the  probability  that  system  j  fails,  given  that  the  aircraft  has  a  malfunction. 
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Once  the  queue  lengths  of  all  network  stations  when  n  customers  in  the  network 
are  determined,  the  next  iteration  of  the  heuristic  can  begin.  The  response  time 
for  all  stations  is  updated  using  the  MVA  algorithm.  Then,  the  cycle  time  for 
each  subnetwork  followed  by  the  system’s  cycle  time  is  calculated.  The  network 
throughput  is  determined  and  finally  utilization  and  queue  length  for  each  station 
is  found  using  the  appropriate  equations,  thus  completing  the  iteration  with  n  -f  1 
customers  in  the  network.  As  before,  this  iterative  process  begins  with  n  =  1  and 
ends  when  n  =  N. 

3,5  Computer  Implementation 

The  MV.'\  heuristic  for  aircraft  sortie  generation  was  programmed  in  Borland’s 
Turbo  Peiscal  for  Windows,  version  1.5.  The  program  calculates  and  displays  the 
mean  values  of  throughput,  queue  length,  response  time  and  server  utilization  for 
each  station  in  the  network.  The  mean  values  of  queue  lengths,  response  time  and 
utilization  at  each  station  are  displayed  for  each  iteration.  The  source  code  for  the 
heuristic,  applied  to  a  network  with  one  fork-join  queue  at  station  4,  is  described 
and  listed  in  Appendix  B.  Appendix  C  contains  sample  output  for  the  program. 


3-15 


IV.  Methodology  Results  and  Comparison  to  LOOM 

4.1  Introduction 

This  chapter  details  the  results  of  the  MVA  heuristic  applied  to  the  bcisic  QNM, 
illustrated  in  figure  1.1  with  figure  1.2  representing  the  concurrent  maintenance  ac¬ 
tivities  of  station  4.  Transition  probabilities  between  queues  are  shown  in  the  figures. 
The  service  time  and  number  of  servers  for  each  queue  are  listed  in  tables  4.1  and  4.2 
below.  The  infinite  number  of  servers  for  nonmaintenance  stations  1  and  2  indicates 
these  are  delay  stations. 


Station 

m 

2 

Fly 

3 

De-Arm 

5 

Turn 

6 

Re- Arm 

Mean  Sve  Time  (hrs) 

0.25 

2.00 

1.25 

1.15 

1.25 

No.  Servers 

00 

00 

2 

8 

9 

Table  4.1,  Service  Time  and  Number  of  Servers  for  Non-Maintenance  Stations 


Fork -Join 

Station 

1 

2 

3 

■ 

5 

6 

Mean  Sve  Time  (hrs) 

2.00 

2.50 

3.25 

3.50 

4.25 

4.50 

No.  Servers 

3 

4 

4 

3 

5 

5 

Table  4.2.  Service  Time  and  Number  of  Servers  for  Fork-Join  Queue  Stations 

After  solving  the  basic  QNM,  specific  network  parameters  were  varied  to  de¬ 
termine  the  sensitivity  of  the  heuristic  to  increased  server  utilization.  At  each  step, 
an  LCOM  simulation  of  the  network  was  created,  the  output  of  which  was  used  as 
a  bciseline  for  determining  the  accuracy  of  the  heuristi' . 
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4-8  LOOM  Simulation  of  the  QNM 


Due  to  the  Air  Force’s  a(.<._ptance  of  the  Logistics  Composite  Model’s  accuracy, 
LCOM  was  used  to  provide  the  baseline  for  determining  the  accuracy  of  the  MVA 
heuristic  for  analyzing  aircraft  sortie  generation.  The  results  of  the  LCOM  simulation 
were  taken  from  the  output  of  the  Matrix  Post  Processor  b  ecause  it  is  the  only 
source  of  LCOM  data  that  reports  aircraft  queueing  time  at  each  station.  Using  this 
information,  plus  the  number  of  aircraft  served  and  the  manhours  worked  at  each 
station,  network  performance  characteristics  for  the  simulation  were  determined. 

The  basic  QNM  was  coded  in  LCOM  version  93. C  format  (9)  and  is  contained 
in  Appendix  A.  As  the  QNM  v/as  modified  for  each  of  the  sensitivity  studies,  the 
LCOM  coded  version  of  the  network  was  correspondingly  changed.  Five  iterations  of 
LCOM  were  run  for  each  network  in  order  to  provide  confidence  limits  on  the  output. 
Of  the  63  simulation  output  variances  calculated  fcr  the  QNM  with  24  aircra.^’t,  only 
six  were  larger  than  0.05,  the  largest  being  0.12.  The  largest  confidence  interval,  as 
a  percentage  of  the  LCOM  output  value,  corresponded  to  the  largest  variance.  This 
occurred  for  the  value  of  the  response  time  for  fork-join  path  5  and  yielded  a  95 
percent  confidence  interval  of  4,28  ±  0.33  hours.  Each  run  of  the  LCOM  simulation 
consisted  of  a  720  hour  warmup  period  followed  by  720  hours  of  data  collection. 

Results  for  the  basic  QNM  with  24  Aircraft 

Tables  4.3  and  4.4  below  compare  the  results  of  the  MVA  heuristic  to  the  output 
of  the  LCOM  simulation  for  the  basic  QNM  with  24  aircraft.  For  all  comparisons 
to  follow,  the  top  line  for  each  station  represents  the  result  of  the  MVA  heuristic 
and  the  bottom  line  represents  the  averaged  output  of  five  LCOM  runs.  The  MVA 
heuristic  applied  to  the  non-maintenance  stations  outside  the  fork-join  queue  yielded 
close  agreement  with  the  LCOM  simulation  output.  Stations  1,  2,  5  and  6.  exhibit 
differences  across  all  performance  characteristics  that  range  only  from  1  to  3  percent. 
Station  3.  the  de-arm  weapons  activity  required  due  to  ground  abort,  had  differences 
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of  approximately  10  percent  between  the  heuristic  results  and  simulation  output  for 
throughput,  queue  length,  and  server  utilization.  However,  the  low-  values  of  the 
performance  characteristics  for  this  station  greatly  reduce  the  significance  of  the 
larger  percentage  difference. 


Station 

Number 

Num  of 
Servers 

Mean 
Svc  Time 

Thru 

Put 

Queue 

Length 

Response 

Time 

Server 

Utilization 

1 

oo 

1.0000 

1.0000 

3.5432 

3.4850 

0.8858 

0.8647 

iilwa 

B 

2 

oo 

0.9500 

0.9528 

3.3661 

3.3203 

6.7322 

6.5311 

6.7322 

6.5311 

3 

2 

gni 

0.1772 

0.1628 

0.2239 

0.2044 

1.2637 

1.2558 

B 

5 

8 

1.1500 

1.1374 

■Vim 

3.5432 

3.4856 

4.1019 

3.9894 

1.1577 

1.1445 

4.0747 

3.9644 

6 

9 

mxiiiM 

3.5432 

3.4844 

4.4426 

4.3608 

B 

B 

Table  4.3.  Comparison  of  Heuristic  and  LCOM  Results  for  Non-Maintenance 
Stations 


While  the  heuristic  result  was  generally  within  5  percent  of  the  simulation  out¬ 
put  for  the  maintenance  stations  on  the  fork-join  queue  paths,  these  differences  were 
less  consistent,  queue  to  queue,  than  for  the  non-maintenance  stations.  For  exam¬ 
ple,  fork-join  queue  path  1  had  only  a  3  percent  difference  in  throughput  between 
the  heuristic  and  simulation,  but  had  a  5  percent  difference  in  response  time  and 
a  10  percent  difference  for  queue  length  and  server  utilization.  On  the  other  hand, 
fork-join  queue  path  2  had  differences  less  than  1  percent  across  all  performance 
characteristics. 


4-4  Sensitivity  Studies 

The  purpose  of  a  sensitivity  study  is  to  determine  the  effect  of  changes  to 
the  QNM’s  parameters  (number  of  customers,  number  of  servers,  service  time)  on 


Station 

Number 

Num  of 
Servers 

Mean 
Svc  I'ime 

m 

Thru 

Put 

Queue 

Length 

Response 

Time 

Server 

Utilization 

1 

3 

— ■ 

0.1742 

0.1686 

ijij  ly 

1.3159 

1.1864 

BUB 

UiilUi 

■lElU 

2 

4 

■wojiiin 

■ 

EBI 

1.8838 

1.8900 

2.6029 

2.5932 

1.8093 

1.8197 

3 

4 

■ 

1.5639 

1.5072 

3.3229 

3.1886 

B 

4 

3 

B 

■ 

1.6862 

1.7386 

B 

5 

5 

B 

0.3500 

0.3525 

m 

4.2644 

4.2843 

1,4877 

1.5056 

6 

5 

4.5000 

4.5466 

0.1435 

0.1478 

ijtey 

ItiSBl 

2.3384 

2.4078 

4.5986 

4.2843 

2.2882 

2.3428 

Table  4.4.  Comparison  of  Heuristic  and  LCOM  Results  for  Fork-Join  Queue 


the  overall  network  throughput.  Any  changes  to  a  specific  station’s  parameters  are 
reflected  in  that  station’s  server  utilization.  A  decrease  in  the  number  of  servers 
or  an  increase  in  service  time  at  a  station  causes  its  server  utilization  to  increase. 
Similarly,  an  increase  in  the  number  of  servers  or  a  decrease  in  service  time  at  any 
station  results  in  lower  server  utilization  at  that  station.  Thus,  any  sensitity  study 
to  changes  in  the  QNM’s  parameters  is  equivalent  to  a  sensitivity  study  to  various 
levels  of  server  utilization. 

In  order  to  determine  the  sensitivity  of  the  MVA  heuristic  to  increased  server 
utilization,  the  basic  QNM  with  24  aircraft  was  re-solved  with  each  of  the  following 
changes  made; 


1.  The  number  of  aircraft  in  the  system  was  doubled  from  24  to  48. 

2.  The  number  of  servers  on  each  of  the  fork-join  queue  paths  Wcls  decreased  to 
one. 

3.  The  number  of  servers  on  fork-join  path  6  was  incrementally  reduced. 
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4-4-^  4^  Aircraft.  The  number  of  aircraft  in  the  system  was  doubled  in  order 
to  increase  overall  congestion,  i.e.,  queueing,  at  each  station.  .Ml  other  parameters 
w’ere  left  unchanged.  As  seen  in  tables  4.5  and  4.6,  while  queue  lengths  and  server 
utilizations  nearly  doubled,  the  results  of  the  heuristic  continue  to  be  in  close  agree¬ 
ment  with  the  output  of  the  LCOM  simulation  for  the  non-maintenance  stations 
outside  the  fork -join  queue.  Differences  between  simulation  and  heuristic  for  these 
stations,  including  station  3,  are  less  than  3  percent  for  all  performance  charac¬ 
teristics.  However,  comparison  of  the  results  for  the  multiple-server  queues  on  the 
fork-join  queue  paths  reveals  differences  that  are  again  inconsistent  in  size  across  the 
queues  and,  for  this  scenario,  are  much  larger.  The  heuristic  result  for  throughput 
at  fork-join  queue  path  1  was  nearly  identical  to  the  simulation  output,  however, 
results  for  queue  length,  response  time  and  server  utilization  differed  by  approxi¬ 
mately  19  percent.  Fork-join  queue  path  4  had  very  small  differences  for  throughput 
and  server  utilization,  but  had  differences  of  approximately  17  percent  for  queue 
length  and  response  time.  These  differences  are  not  attributable  to  the  variance  of 
the  LCOM  simulation  output;  each  of  the  large  errors  noted  above  lies  well  outside 
the  95  percent  confidence  interval  for  the  corresponding  LCOM  simulation  output. 

We  note  that  the  largest  errors  occurred  on  the  two  fork-join  queue  paths  that 
have  server  utilization  levels  above  80  percent  and  that,  ai  levels  of  server  utilization 
75  percent  and  below,  the  MVA  heuristic  results  were  within  a  few  percentage  points 
of  the  LCOM  simulation.  This  leads  us  to  believe  that  the  MVA  heuristic  loses 
accurracy  for  stations  on  the  fork-join  queue  paths  when  server  utilizations  are  high. 

4-4-^  Single-Server  Queues  on  Fork-Join  Queue  Paths.  In  an  attempt  to 
characterize  the  levels  of  server  utilization  at  which  the  heuristic  no  longer  pro¬ 
vides  accurate  approximations  of  network  performance  characteristics,  a  pathologi¬ 
cal  maintenance  scenario  Wcis  created.  Using  the  basic  QNM  with  24  aircraft,  the 
number  of  servers  on  each  of  the  fork-join  queue  paths  was  reduced  to  one.  That  is, 
each  of  the  multiple-server  queues  was  replaced  by  a  single-server  queue.  Also,  in 
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Table  4.5.  Comparison  for  Non-Maintenance  Stations,  48  Aircraft 


Station 

Number 

Num  of 
Servers 

Mean 

Sve  Time 

■ 

Thru 

Put 

Queue 

Length 

Response 

Time 

Server 

Utilization 

1 

3 

2.0000 

1.9634 

0.1742 

0.1673 

1.0292 

1.0253 

2.9612 

2.4781 

2.8773 

2.4161 

2.0583 

2.4781 

2 

4 

2.5000 

2.4946 

Ml 

Infill 

m 

3.5643 

3.2673 

2.5502 

3.0886 

3 

4 

Wiil 

3.9132 

3.8514 

4 

3 

■ 

5.1510 

4.2189 

2.5118 

2.5444 

5 

5 

4.2500 

4.2789 

IWal 

0.5836 

0.6011 

Mifga 

4.4402 

4.5148 

2.4802 

2.5692 

6 

5 

4.5000 

4.3628 

■ 

||Ky^ 

5.0409 

4.9747 

3.8150 

3.8906 

Table  4.6.  Comparison  for  Fork-Join  Queue,  48  .Aircraft 
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order  to  keep  the  confidence  interval  on  the  simulation  output  reasonably  tight,  20 
iterations  of  the  LOOM  simulation  were  conducted. 

Table  4.7  indicates  that  the  heuristic  continued  to  provide  highly  accurate  re¬ 
sults,  within  1  to  2  percent  of  the  LCOM  simulation,  for  non-maintenance  stations. 
As  shown  in  table  4.8,  the  heuristic  also  accurately  calculated  throughput  and  server 
utilization  for  the  fork-join  queue  stations.  However,  for  the  performance  character¬ 
istics  of  queue  length  and  response  time,  the  differences  in  results  increeised  signifi¬ 
cantly.  The  largest  differences  did  not  correspond  to  the  highest  server  utilization, 
the  single-server  queue  on  fork-join  queue  path  6  with  98  percent  server  utilization. 
Instead,  the  heuristic  was  most  inaccurate  for  fork-join  queue  path  4  where  results 
differed  from  the  LCOM  simulation  by  over  50  percent  while  utilization  was  only  59 
percent.  Fork-join  queue  paths  1,  2,  3  and  5  had  differences  between  the  heuristic 
and  simulation  that  ranged  from  21  to  33  percent  and  server  utilizations  that  ranged 
from  52  to  78  percent. 

The  apparent  accuracy  of  the  heuristic’s  results  for  fork-join  path  6  is  surprising 
for  two  reasons.  First,  the  heuristic’s  results  for  queue  length  and  response  time  are 
below  the  values  of  the  LCOM  simulation  output,  which  have  95  percent  confidence 
of  11.25  ±  0.78  aircraft  and  51.95  ±  4.43  hours,  respectively.  For  the  stations 
on  fork-join  queue  paths  1  through  5  the  heuristic’s  results  for  queue  length  and 
response  time  are  well  above  their  corresponding  LCOM  output  values.  Second,  the 
server  utilization  on  fork-join  queue  path  6  is  98  percent,  by  far  the  highest  server 
utilization,  while  the  heuristic’s  results  for  the  path  are  the  most  accurate.  The  final 
sensitivity  study  will  explore  the  accuracy  of  the  MVA  heuristic  when  the  server 
utilization  rate  is  increased  for  only  one  fork-join  queue  path,  path  6  in  particular. 

4-4-3  Incremental  Reduction  of  Servers  on  Fork-Join  Path  6.  In  order  to 
study  the  effects  of  high  server  utilization  on  a  single  fork-join  queue  path,  a  final 
sensitivity  study  was  conducted  on  the  48  aircraft  scenario.  The  number  of  servers 
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Table  4.7.  Comparison  for  Non- Maintenance  Stations,  One  Server  per  Fork-Join 
Path 


Station 

Number 

Num  of 
Servers 

Mean 
Svc  Time 

m 

Thru 

Put 

Queue 

Length 

Response 

Time 

Server 

Utilization 

1 

1 

2.0000 

1.9949 

DBEI 

0.2633 

0.2636 

11 

4.1791 

3.4596 

0.5266 

0.5246 

2 

1 

2.5000 

2.4896 

■ 

0.3087 

0.3164 

3.1037 

2.5591 

10.0537 

7.8680 

IHi 

3 

1 

3.2500 

3.1482 

H 

0.2008 

0.2074 

n 

4 

1 

0.1215 

0.1164 

0.1836 

0.1792 

BBS 

9.5397 

6.3453 

0.6427 

0.5924 

5 

1 

0.0988 

0.1021 

0.1493 

0.1569 

1.6943 

1.4312 

11.3471 

9.0352 

■■ 

6 

1 

i 

0.1435 

0.1429 

0.2169 

0.2194 

10.0126 

11.2474 

46.1607 

51.9500 

■I 

Table  4.8.  Comparison  for  Fork-Join  Queue,  One  Server  per  Path 
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on  fork-join  queue  path  6  was  incrementally  reduced  by  one  thereby  incrementally 
increasing  server  utilization  at  this  station.  All  other  parameters  in  the  network  were 
left  unchanged.  At  every  station  in  the  QNM  other  than  the  station  on  fork-join 
queue  path  6,  the  MVA  heuristic  accurately  determined  the  network  performance 
characteristics,  differing  at  most  by  5  percent  from  the  LCOM  sim>ilation  output. 
Table  4.9  displays  for  fork-join  queue  path  6  the  results  of  the  MVA  heuristic  and 
the  averaged  output  of  5  LCOM  runs  for  each  level  of  servers. 


Num  of 
Servers 

Thru 

Put 

Queue 

Length 

Response 

Time 

Server 

Utilization 

5 

lifMI 

5.9461 

5.5705 

3.8150 

3.8906 

4 

0.8132 

0.8322 

8.2908 

8.4686 

10.2068 

10.1976 

3.6553 

3.8622 

3 

0.6586 

0.6797 

■Mjllgl 

25.8409 

26.2464 

2.9637 

2.9753 

2 

25.2010 

28.6392 

57.1198 

65.9732 

IH 

1 

H 

32.3064 

34.3025 

146.2026 

163.6027 

Table  4.9.  Comparison  for  Fork-Join  Queue  Path  6  Study 


MVA  heuristic  accurately  determines  the  throughput  and  server  utilization  as 
the  number  of  servers  decreased.  The  heuristic’s  result  for  response  time  is  7  percent 
above  the  LCOM  simulation  output  when  there  are  5  servers  on  the  path  and  server 
utilization  is  at  78  percent.  As  the  number  of  servers  decreases  and  server  utilization 
increases  above  95  percent,  the  heuristic’s  calculation  for  the  path’s  response  time 
moves  to  13  percent  below  the  LCOM  simulation  output.  The  analytical  value  for 
queue  length  follows  a  similar  trending  as  server  utilization  increases,  reaching  a  max¬ 
imum  difference  of  approximately  10  percent.  These  observations  were  duplicated 
for  fork-join  queue  path  4,  including  the  incident  of  high  accuracy  at  approximately 


95  percent  utilization.  The  reasons  for  this  behavior  at  the  very  highest  levels  of 
server  utilization  remain  unclear. 
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V.  Conclusions  and  Recommendations  for  Future  Research 

5.1  General  Observations 

This  effort  was  motivated  by  the  Air  Force’s  current  policy  of  using  LCOM  to 
determine  maintenance  manpower  requirements.  Preparing,  running,  and  optimiz¬ 
ing  a  new  Logistics  Composite  Model  simulation  is  a  major  undertaking  requiring 
months  of  effort.  System  failure  data  must  be  collected,  detailed  task  structures 
must  be  determined  and  validated,  aircraft  mission  scenarios  and  equipment  config¬ 
urations  must  be  identified,  LCOM  code  must  be  written  and  debugged  and,  finally, 
simulation  output  must  be  analyzed.  The  manpower  study  itself  requires  multiple 
runs  of  the  simulation  in  order  to  arrive  at  the  optimum  mix  of  resources  to  support 
the  desired  mission  scenario  (29). 

While  the  methodology  developed  in  this  study  is  a  lumristic,  it  quickly  pro¬ 
vides  a  reasonably  accurate  estimate  of  maintenance  manpower  requirements  for  a 
given  aircraft  sortie  generation  process.  The  data  required  to  build  the  QNM  for  the 
heuristic  is  less  complicated  than  the  data  used  in  LCOM.  Individual  system  failure 
data  is  not  needed;  analysts  must  only  determine  the  probability  that  each  mainte¬ 
nance  speciality  works  on  an  aircraft  that  has  some  malfunction.  The  average  crew 
size  and  job  completion  time  for  each  maintenance  speciality  replaces  detailed  ttisk 
structures.  All  of  this  data  can  be  eaisily  determined  from  the  information  collected 
by  the  Air  Force  Job  Control  function  that  monitors  maintenance  activities  for  the 
actual  aircraft  sortie  generation  process  to  be  studied. 

The  comparisons  of  the  results  of  the  MVA  heuristic  to  the  outputs  of  the 
LCOM  simulation  indicate  that  the  heuristic’s  accuracy  decreeises  as  server  utiliza¬ 
tion  increases.  When  server  utilization  is  kept  in  realistic  ranges  (utilization  of 
Air  Force  maintenance  personnel  rarely  exceeds  70  percent  (29)),  the  results  of  the 
heuristic  for  a  single  fork-join  queue  are  very  accurate.  For  non-maintenance  sta¬ 
tions,  results  are  within  1  to  2  percent  of  the  LCOM  simulation  output.  For  stations 
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on  the  fork-join  queue  paths,  heuristic  results  are  within  5  percent  of  the  LCOM 
simulation’s  output  for  that  portion  of  the  network. 

5.8  Limitations  of  the  Methodology 

A  key  step  in  the  M\"A  heurisitc  for  analysis  of  aircraft  sortie  generation  is  the 
averaging  of  all  of  the  possible  network  cycle  times  ^  o  get  the  actual  network  cycle 
time.  For  our  QNM  with  only  one  fork-join  queue  that  has  6  paths,  this  entailed  the 
consideration  of  only  7  cases.  However,  if  we  desire  to  expand  this  methodology  to 
complicated  QNMs  that  have  more  than  one  fork-join  queue,  this  averaging  process 
could  quickly  become  computationally  tedious.  For  example,  a  QNM  with  three 
fork-join  queues  with  3,4  and  6  paths  would  entail  the  calculation  and  averaging  of 
140  different  cycle  times. 

5.3  Recommendations  for  Future  Research 

The  M\^A  heuristic  developed  in  this  study  is  a  first  step  towards  the  creation 
of  an  analytical  model  that  accurately  relfects  a  system  that  has  heretofore  only  been 
modeled  via  simulation.  Additional  work  is  required  to  increase  the  model’s  fidelity. 
The  topics  suggested  below  are  further  steps  in  that  direction. 

5.3.1  Heuristic  Effect  of  Assumption  (Al).  The  methodology  developed  by 
Rao  and  Suri  (26)  and  extended  here  hinges  on  the  hueristic  application  of  the 
Arrival  Theorem  to  fork-join  queues.  Assumption  (Al)  is  exact  when  there  is  only 
one  customer  in  the  network  and  is  approximate  for  more  than  one  customer.  The 
heuristic  effect  of  the  assumption  and  its  relation  to  network  parameters  such  as  the 
number  of  customers  and  the  rumber  of  fork-join  queue  paths  should  be  investigated. 
Ideally,  this  investigation  coui  'i  if:ad  to  an  adjustment  scheme  that  could  be  applied 
to  network  performance  characteristics  in  order  to  increase  solution  accuracy. 


5.3.2  Heuristic  Effect  of  Assumption  (A2).  The  work  presented  in  this  study 
is  also  dependent  upon  Rao  and  Suri’s  second  assumption,  that  the  response  time  of 
an  item  at  fork-join  queue  path  i  is  exponential.  This  assumption,  I'ke  Assumption 
(Al),  is  exact  only  when  there  is  only  one  customer  in  the  network  and  is  approx 
imate  for  more  than  one  customer.  Making  this  assumption  allows  us  to  apply  the 
memoryless  property  of  the  exponential  distribution  in  order  to  determine  the  length 
of  time  a  sibling  customer  on  a  fork-join  queue  path  must  wait  for  its  other  siblings 
before  being  re-joined.  The  actual  distribution  of  the  r  sponse  time  is  gamma  with 
parameters  dependent  upon  the  number  of  customers  in  the  queue.  The  gamma 
distribution  hcis  a  smaller  variance  than  the  exponential,  thus  the  heuristic  effect  of 
the  assumption  is  to  overestimate  the  response  time.  The  size  of  this  effect  and  its 
relation  to  network  parameters  such  a^  the  number  of  customers  and  the  number  of 
fork-join  queue  paths  should  ■  /estigated.  This  investigation  could  also  lead  to 
an  adjustment  scheme  that  cou)d  be  applied  to  network  performance  characteristics 
in  order  to  increase  solution  accuracy. 

5.3.3  Machine' Operator  Interference.  Several  simplifying  assumptions  (de¬ 
tailed  in  Chapter  1)  were  made  in  this  study  that,  while  making  the  QNM  easier 
to  model  and  solve,  reduced  model  accuracy.  One  simplification  weis  the  averaging 
of  crew  sizes  and  service  times  over  each  maintenance  speciality’s  tasks.  Another 
dealt  with  ihe  machine-operator  interference  caused  by  separate  queues  using  the 
same  resource.  These  two  issues  can  be  directly  addressed  without  simplification  by 
iHodeling  the  aircraft  sortie  generation  process  with  a  hierarchal  system  of  QNMs 
that  exhange  information  (28:267). 

5.3.4  Validation  against  an  Actual  Scenario.  In  oder  to  assess  the  mathe¬ 
matical  accuracy  of  the  heuristic,  the  validation  presented  in  this  study  compared 
the  heuristic  results  to  the  output  of  an  LOOM  simulation  of  the  exact  same  queu^ 
ing  network  model.  However,  the  suitability  of  using  the  heuristic  as  a  complement 
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to  LCOM  was  not  addressed.  A  validation  for  this  purpose  would  compare  an  ex¬ 
isting  LCOM  simulation  of  an  actual  aircraft  sortie  generation  process,  of  which 
many  suitable  examples  exist,  to  the  results  of  the  heuristic  applied  to  an  equivalent 
QNM.  Due  to  the  complexity  of  LCOM  tasks  structures,  this  “real-world”  valida¬ 
tion  should  be  performed  as  a  subset  the  machine-operator  interference  research 
proposed  above. 


Appendix  A.  Sample  LOOM  Inpxd 
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Appendix  B.  Program  MVAHEUR 


B.  1  Description 

The  MVA  heuristic  for  analysis  of  aircraft  sortie  generation  is  implemented 
in  the  program  MVAHEUR.  The  program  is  written  in  Borland’s  Turbo  Pascal  for 
Windows,  version  1.5  and  is  easily  changed  into  standard  Pascal  format  so  that  it 
can  be  run  on  any  mainframe  or  personal  computer  with  a  Paiscal  compiler.  The  pro¬ 
gram,  listed  below,  contains  substantial  internal  documentation  including  variable 
definitions  and  brief  descriptions  of  each  subroutine.  An  overview  of  the  relationship 
between  the  various  program  components  is  illustrated  in  figure  B.l. 

B.2  Data  Files 

To  use  program  MVAHEUR  data  should  first  be  entered  into  a  text  file  called 
MVA. DAT  in  accordance  with  the  format  displayed  in  Table  B.2.  The  general  format 
is  paralleled  by  the  data  file  for  the  QNM  example  with  24  aircraft  solved  in  Chap¬ 
ter  4.  In  the  data  file,  N  is  the  number  of  aircraft  and  M  is  the  number  of  stations 
(including  the  fork-join  queue  as  station  4  but  not  including  the  stations  on  the 
fork-join  paths).  ASeruer[i]  is  the  number  of  servers  at  each  of  these  stations  with 
NServer[A]  being  any  value.  5[z]  is  the  service  time  at  each  of  the  stations  with, 
again,  S[4]  being  any  value.  The  matrix  of  transition  probabilities,  P[i,  j],  is  followed 
by  the  fork-join  queue  parameters  for  each  path.  These  are  the  number  of  servers 
on  each  path:  F J Server[i],  the  service  time:  FJ5[z],  and  the  path  probabilities: 
FJP[i].  Output  of  the  program  is  written  to  a  text  file  called  MVA. OUT.  The  output 
file  for  the  data  file  below  (the  QNM  with  24  aircraft)  is  contained  in  Appendix  C. 
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Figure  B.l.  Relationship  Between  Components  of  Program  MVAHEUR 
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Table  B.l.  Data  Format  for  Program  MVAHEUR 
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B.S  Program  Listing 


program  NVAHBUR; 


*> 

Appllas  a  Naan  Value  Analyais  Huariatic  to  datarmina  a> 

{a  parforaanca  maasuraa  for  atationa  in  a  cloaad  network  *} 
{a  of  aultiaarvar  quauaa  where  one  queue  ia  a  Fork-Join  *} 
{a  Queue  a} 

{a  a> 

{a  Language:  Turbo-Pascal  for  Windows  l.S  a} 

{a  •} 

{a  Source:  R  C  Jenkins,  rjenkinsflsifit.af .ail.  Box  4498  a} 

{a  AFIT/EBS,  Wright-Pattaraon  AFB  OH  46433-7765  a} 

{a  *y 

i**m***********^**************»»***********t*******m*^*****} 
uaea  VinCrt; 


conat 

Imax  =  48;  {maximum  number  of  aircraft  in  network} 
Mmax  =  9;  {maximum  number  of  stations  in  network} 
Paths  =  6;  {maximum  number  of  paths  in  the  FJ  queue} 

type  Age  =  (old, new); 

NIntArray  =  array  Cl ■. Mmax]  of  integer; 

NRealArray  »  exray  Cl.  .Mmax]  of  real; 

PNatrix  arrayCl .  .Mmax,  1.  .Mmax+l]  of  real; 

FJIntArray  s  array  Cl. .Paths]  of  integer; 

FJRealArray  =  array  Cl .  .Paths]  of  real; 

OutNatrix  =  arrayCl.  .Nfflax,0.  .Vmax]  of  real; 
FJOutNatrix  =  arrayCl .  .Patha,0.  .Imax]  of  real; 
■RealArray  =  arrayCl.  Vmax]  of  real; 

FJTimeArray  =  axT...yCC.  .Path”]  of  real; 

FJTlmeArrayS  -  array CO. .Paths^l]  of  real; 

SubPrArray  =  arrayCl . .Paths+l]  of  real; 

CTMatrix  =  arrayCl .  ■Paths'*-1,0.  .Mmax]  of  real; 


var  I,J,K,FJ;  integer;  {counters} 

M,l:  integer;  {number  of  stations,  customers  in  QHH} 

CycleTlme:  real;  {cycletime  of  QHH} 

Lambda:  IRealArray;  {throughput  of  QIN} 

IServer:  NIntArray;  {number  of  servers  at  each  station} 

FJServer:  FJIntArray;  {number  of  servers  at  fork-join  queue  path} 


{service  time  of  stations} 

{visit  ratios  of  stations} 

{service  time  of  fork-join  queue  path  stations} 
{fork-join  queue  path  probabilities} 

{transition  probability  matrix  for  QNM} 

C:  arrayCl ■ .Nmax,l. .Imax]  of  integer;  {number  of  busy  servers  per  station} 
FJC:  array  Cl . .Paths , 1 . .Imax]  of  real;  {number  of  busy  servers  on  fj  paths} 


S:  NRealArray; 

V:  NRealArray; 

FJS ;  FJRealArray ; 
FJP:  FJRealArray; 
P:  PMatrix; 
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Q.R.U:  OutM&trix;  {output  matricsB  for  uon-maintenanco  stations} 
FJR,FJQ,FJU:  FJOutNatrix;  {Output  data  matricsa  for  fork-join  queua} 
Pr:  array [old .. non]  of  OutNatriz;  {probabilitlos  for  marginal  local...} 
FJPr;  array [old. .non] of  FJOutNatrix;  {...bedanco  thoorem  application} 
Dat.Out:  tozt;  {input  and  output  filos} 


procoduro  Error(ErrCode,I,K:  intogor); 

{Reports  orrors  -  callod  by  program  MVAHEUR  and  procoduro  RoadData} 

bogin 

nritoln; 
caso  ErrCodo  of 

1:  nritolnC' ERROR:  H  not  in  {1 . . * ,Hmax:2, *}' ) : 

2:  nr itolnC 'ERROR:  K  not  in  {l. . ',Naax:2, ; 

3:  nritolnC 'ERROR:  routing  probs  from  station  ’,1:2, 

'  do  not  sum  to  1.0000'}; 

4:  nritolnC'WARHIIG:  Pr{' ,1:2, ' ,0, ’ ,K:2, ’}<0;  ', 

'may  have  numorical  probloms'); 

6:  nritolnC'FJWARIIIG:  Pr{' .1:2, ' ,0. ' ,K;2, '}<0;  ', 

'may  have  numorical  problems’); 
end;  {case} 

if  (ErrCodoM  or  5)  then  bogin 

nritolnCPr'-as  <ontor>  to  resume’); 
roadln ; 

end  else  begin 

nritolnC' Program  terminated;  press  <onter>  to  exit’); 

readln; 

halt; 

end;  {else} 
end;  {Error} 
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-  procedure  ReadData(var  IiN:  integer; 

var  HServer:  Mint Array; 
var  FJServer:  FJInt Array; 
var  S:  MRealArray; 
var  FJS,FJP:FJRealArray; 
var  P;  PNatriz; 
var  Dat;  text); 

{r.eads  input  data  from  lile  assigned  to  text  variable  'Dat’} 

{Called  by  program  MVAHEUR} 

var  I,J:  integer;  {counters} 

Psum:  real;  {holds  sum  of  transition  probabilities  matrix  rows} 
begin 

readln(Dat,l,N) ; 

if  not  (I  in  Cl..Imax])  then  Error(1.0,0) ; 
if  not  (M  in  [l..Mmaz])  then  Error(2.0.0) ; 
for  I:=l  to  M  do  read(Dat,IServerCl}); 
for  I:=l  to  M  do  read(Dat .SCI] ) ; 
for  I:=l  to  H  do  begin 
Psua;30; 

for  J:3l  to  N  do  begin 
read(Dat,P[J,I]); 

Poum:*Poum+PCJ,I] ; 
end;  {for} 

if  (Psum<0.9999)  or  (Psum>1.0001)  then  ErrorO.I.O) ; 
end;  {for} 

for  1:=  1  to  Paths  do  read(Dat,FJServer[I]} ; 
for  I:=  1  to  Paths  do  read(Dat,FJSCl]) ; 
for  I:=  1  to  Paths  do  read(Dat,FJP[I]); 

P[l,l]:=l; 

for  J:=2  to  M  do  P[1,J]:=0; 

P[1,M+1];=1; 
for  I:  =2  to  M  do  begin 
PCI. I]  :=PCI.I]-1; 

PCI.M+1]:=0; 
end;  {for} 
end;  {ReadOata} 
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proc«dar«  VSolv«(N;  integer: 

var  P:  PKatrix; 
var  V:  NRealArray); 

{Solves  simultaneous  equations  to  obtain  visit  ratios} 
{Called  by  KVAHEUR} 

var  I,J,K,IC,KK,NN,IT.IS:  integer;  {counters} 

B,V,C:  real; 

ID:  Mint Array; 

Y:  MRealArray; 

begin 

MN:=N-^1; 

for  I:=l  to  M  do  IDCI] :=!; 

K;=l; 

repeat 

KK;sK+l: 

IS:=>K; 

IT:=K; 

B:=abs(PCK.K]): 

for  I:=K  to  M  do  for  J:=K  to  M  do 
if  (absCPCi , J] )>B)  then  begin 
IS:=I; 

IT:=J; 

B:=ab8(PCl,J]); 
end;  {if} 
if  (IS>K)  then 

for  J ; =K  to  NN  do  begin 
C:=P[IS,J]; 

P[IS,J]  :=P[K,J]; 

PCK,J]:=C; 
end;  {for} 

if  (IT>K)  then  begin 
IC:=IDCK]  ; 

IDCK]  ;=IDCIT]  ; 

ID [IT] :=IC; 

for  I:=l  to  K  do  begin 
C;=P[I.IT]: 

PCI.IT]:=P[I,K]; 

PCI.K] :=C; 
end;  {for} 
end;  {if} 

for  J;=KK  to  MM  do  begin 
P[K.J]:=PCK,J]/P[K.K]; 
for  I:=KK  to  N  do  begin 
W:=P[I,K]*PLK,J]; 

PCI.J] :=PCl,J]-tf: 
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il  (abs(P[I,J])<0.00001*ab8(W))  than  P[I. J] :=0 
and;  {{or} 
end;  {lor} 

K:=KK: 

until  (K=M); 

Y[M] :=P[M,HM]/P[M,M]; 
lor  I:=l  to  M-1  do  begin 
K;=M-I: 

KK:=K+1; 

Y[K]:=P[K.MM]  ; 

for  J:=XK  to  M  do  Y [K]  :  =Y [K] -P DC . J] *Y [J]  ; 
and;  {lor} 

lor  I:=l  to  M  do  lor  J:=l  to  M  do 
il  (ID[J]=I)  than  V[I];=Y[J]; 
and;  {VSolva} 
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{m0*****it***************t*************************0t*****t*y 


procedara  AddTemCJ.Bign:  integer; 

nua, den: real; 
var  I:  integer; 
var  pr:  real; 
var  T:  FJRealArray) ; 

{calculates  terns  lor  fork-join  queue  waiting  probabiities} 
{Called  by  procedure  FJPTiae  and  calls  itself} 

var  K;  integer;  {counter} 

begin 

num:=nun+TCJ] ; 
den:=den+T[J]  ; 

8ign:=8ign*(-l) ; 
pr ;  =pr'i-8  ign*nun/den ; 
if  (J<I-1)  then 

for  K:=J+1  to  I-l  do  AddTerm(K, sign. nun, dan, 1, pr, T) ; 
end;  {AddTem} 

{eeeeeeeeee**ee«****ee«***e«****e«eeee*«*ee*««eee*e**ee*ee*} 

procedure  PJPTine  (  Paths:  Integer; 

T:  FJRealArray; 

FJR:  FJOutmatriz; 

FJP:  FJRealArray; 
var  Tine:  Real); 

{detemines  response  time  of  fork-join  queue  for  each  path} 
{Called  by  CalcCycleTime,  calls  AddTem} 

var  FJTPr.Tem:  FJRealArray: 

VTime:  FJTimeArray; 
sign.l.J:  integer; 
nun. den, pr:  real; 

begin 

for  1:=1  to  Paths  do  begin 
pr:=0; 
nun:sO; 
den:=T[I] ; 

lor  J:=l  to  I-l  do  begin 
8ign:=-l; 

AddTora(J  .sign.iium.den.I.pr  ,T) ; 
end;  {lor} 

FJTPrCl]:=pr; 
end;  {lor} 
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lor  I:=l  to  Paths  do  TemCl]  :=FJTPrCl3*FJRCl,K]*FJP[I]  ; 
WTia®CO]:=0: 

lor  I:=  1  to  Paths  do  WTimaCi]  :=WTia6[I-l]'*-TenB[I]  ; 
TiBs:s  FJRCl.K]  *  UTias[6]; 
end;  {FJTine} 
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procoQura  CulcCy ‘'.laTinia (Paths, M:  Integer; 

V:  HRealArray; 

FJR:  FJOutDatrix; 

FJP:  FJRealArray; 

R:  OutHatrlx; 

var  CvcleT'.Be:  Rsal); 

{averages  Icrk-join  response  times  for  all  paths} 

{Called  bv  program  MVAHEUR,  calls  FJPTime} 

var  PTime:  FJRealArray;  {resporse  time  ol  subnetvorks} 

T;  FJRealArray;  {inverse  fork-join  queue  response  times} 
TotProb:  FJTimeArray;  {sums  of  subnetuork  probabilities} 
TotTine:  FJTiroArray^;  {sxuns  of  cycletimes  adjusted  by  probe} 
I,J:  integer;  {counters} 

Time:  Real;  {re8''onBe  time  of  fork- join  queue} 

CT:  CTMatrix;  {cycletime  for  each  subnetwork} 

SubPr;  SubPrArray;  ^probability  of  a  subnetuork} 

begin 

for  I:=  1  to  Paths  do  begin 
TC1]:=1/FJRCI,K3; 
end;  {for} 

FJPTimeCPaths ,T, FJR, FJP, Time) ; 

PTime [l]:=Time: 

for  J:si  to  6  do  begin 

for  1:=!  to  6  do  begin 
TCI];=TCI+1]; 

FJR[I,K] :=FJRr>+i,K]; 

FJP [I] :=FJPCI+1] ; 
end;  {for} 

TC6];=T[1]; 

FJR[6,K] :=FJR[l,K] ; 

FJP [6] :=FJP[1] ; 

FJPTimeCPaths, T, FJR, FJP, Time); 

PTime [J+i] : =Time ; 
end; {for} 

For  I:=l  to  Paths  do  begin 
R[4.K] :=PTime[I] ; 

CT[I,0];=0; 

for  J;=l  to  M  do  CT[I.J] :=CTCl,J-l]+V[J]*R[J,K] ; 
end;  {for} 

CT  [7  ,  M]  :  =  V  [  1  ]  *R  [  1 , K]  +  V  [2]  eR [2 , K]  +  V  [6]  * R [6 . K]  +  V  [6]  ♦  R  [6 , K]  ; 


TotProbCO] :=0; 


lor  I:=l  to  Paths  do  begin 

TotProbCl] :=  TotProbCl-l]+FJP[l3*V[4] ; 
end;  {lor} 

lor  I;=l  to  Paths  do  SubPr[I] :=FJPCI]*V[4] ; 

SubPrC7]:=l-V[4]; 

TotTimeCO] :=0; 

lor  I:=l  to  Paths-t-l  do  begin 

TotTime  [I]  :  =Tot t  ime  [I- 1]  +CT  [I .  M]  eSubPr  [I]  /  (Tot Prob  [6]  +  ( 1  -V  [4]  ) ) ; 
end;  {lor} 

CycleTim«:=TotTime[7] ; 
end;  {CalcFJTime} 
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procedure  WrlteOut(R,M:  integer; 

var  X;  OutKatrix; 
var  Out:  text); 

{Vritee  conplete  output  matrix  lor  perlormance  measures} 
{Called  by  procedure  Report} 

var  I,K:  integer;  {counters} 

begin 

eriteln(Out) ; 
write (Out,'  H'); 

lor  I:  =  l  to  M  do  write(0ut,’  i=’,I:'l); 
writeln(Out) ; 
write (Out,’  — ’); 

lor  I:=l  to  M  do  write(Out,’  - ’); 

writeln(Out) ; 
lor  K;=l  to  H  do  begin 
write(0ut,K:3) ; 

lor  I;=l  to  M  do  write(Out,XCl,K] :7:3) ; 
vriteln(Out) ; 
end ;  {lor} 
wrlteln(Out) ; 
writeln(Out) ; 
end;  {VfriteOut} 


. . . 

procedure  FJVriteOut(H , Paths :  integer: 

var  FJX:  FJOutMatriz; 
var  Out:  text); 

{Writes  output  lor  Fork-Join  Queue  perlormance  measurej} 
{Called  by  procedure  FJReport} 

var  I,K:  integer;  {counters} 

begin 

writelu(Out) ; 
write (Out,’  H’); 

lor  I:=l  to  Paths  do  write(0ut,’  i=’,I;2): 
writeln(Out) ; 
write (Out,'  — '); 

lor  I:=l  to  Paths  do  wiite(0ut,’  - '); 

writeln(Out) ; 
lor  K:=l  to  X  do  begin 
write(0ut,K;3) ; 

lor  I:=l  to  Paths  do  wri*’e(Out,FJX[I,K]  :7:3) ; 


B-13 


writelnCOut) ; 
•nd;  {lor} 
vriteln(Out) ; 
«rlt«ln((hit) ; 
end;  {FJWritePvit} 


procedure  Repoic(H,H:  integer; 

Lambda:  HRealArray; 
var  HServer:  Mint Array; 
var  S,V:  MRealArray; 
var  QjR,U:  OutHatrix; 
var  Out:  text); 

{Writes  output  report  to  lile  assigned  to  text  variable  'Out'> 
{Called  by  MVAHEUR,  calls  procedure  WriteOut} 

var  I:  integer;  {counter} 

SLambda:  real;  {individual  station's  throughput} 

begin 

vriteln(Out, 'NVA  PERFORMAHCE  REPORT’); 

«riteln(Out) ; 

vriteln (Out. 'Humber  of  Customers  (I)  =’,1:3); 
trriteln(Out) ; 

mritelnCOut, 'Stat  Hum  of  Mean  Visit  Thru  ’, 

'  Queue  Respns  Server'); 

«riteln(Out, '  (i)  Servrs  Svc  Tm  Ratio  put  ', 

'  Length  Time  Util  '); 

writelnCOut,  ’ -  -  -  -  - ', 

’  . -  - '): 

for  I;=l  to  N  do  begin 

SLambda : =Lambda [H] eV [1] ; 

wr iteln(0ut , I : 3 , HServer [I] : 7 , S [I] : 10 : 4 , V [I] : 8 : 4 , 
SLafflbda:8:4.qCl,H] :8:4,RCl,H] :8:4,UCl,H] ;8:4) ; 
end;  {for} 

«riteln(Out) ; 

«riteln(Out) ; 

vriteln (Out, 'Average  Queue  Lengths  (including  service)'); 
WriteOut(H,M,Q,Out) ; 

sriteln(Out, 'Average  Response  Times  (including  service)'); 
WriteOut(H,M,R,Out) ; 

Vriteln(Out, 'Average  Utilizations’); 

WriteOut (N,M,U, Out) ; 
end;  {Report} 
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procedure  F J Report (R,N, Paths:  integer; 

V :  MRealArray ; 

Lambda:  MRealArray; 
var  FJServer:  FJlnt Array; 
var  FJS.FJP:  FJRealArray; 
var  FJq.FJR.FJU:  FJOutMatrix; 
var  R:  OutNatrlx; 
var  Out :  text) ; 

{Writes  output  tor  FJ  queue  to  file  assigned  text  variable  'Out'} 
{Called  by  MVABEUR.  calls  FJUriteQut} 

var  I:  integer;  {counter} 

FJLambda:  real;  {individual  fork-join  path  throughput} 
begin 

vriteln(Jut,’MVA  PERFORMABCE  REPORT  FOR  FORK-JOII  QUEUE’); 
Briteln(Out) ; 

sr it eln (Out, 'Humber  of  Customers  (H)  =’,M:3); 
vriteln(Out) ; 

vriteln(Qut, *Stat  Hum  of  Mean  Visit  Thru 

'  Queue  Respns  Server ’ )  ; 

»riteln(Out, ’  (i)  Servrs  Svc  Tm  Ratio  put 
'  Length  Time  Util  ’ ) ; 

«riteln(Out, ' -  — —  -  -  — - 

'  . ’): 

for  1:-1  to  Paths  do  begin 

FJLambda :=Lambda[H]eFJPCl]*V [4] ; 

vriteln(0ut,I:3,FJServerCl] :7,FJSEl] : 10:4,  V[4]*FJP[I] :8:4, 
FJLambda: 8 :4,FJQ[I,H] :8:4.FJR[I,H] :8:4,FJU[I,H] :8:4) ; 
end;  {for} 
vriteln(Out) ; 
vriteln(Out) ; 

«riteln(Out, 'Average  Fork-Join  Queue  Lengths  (including  service)’) 
FJWritaOut(H, Paths, FJQ, Out) ; 

vriteln(Out, 'Average  Reeponse  Times  (including  service)'); 
FJWriteOut(M, Paths, FJR, Out) ; 

Vriteln(Out, 'Average  Utilizations') ; 

F JWr it eOut ( H , Paths , FJU , Out ) ; 
end;  {FJReport} 
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bagin 

assignCDat. 'MVA.DAT'); 
a88ign(0ut, 'MVA.OUT'): 
rasatCDat) ; 
reBTite(Out) ; 

Britain; 

Britaln(’***  Running  program  MVA  ••*'): 

Britain; 

Britaln('>  Reading  data  ...'); 

ReadData( H , H , IServar , F JSarver .S,FJS,PJP,P, Dat ) ; 

Britaln('>  Calculating  visit  ratios 
VSolve(H,P,V); 

vritaln(‘>  Calculating  parfomance  maasuras 
lor  I;=l  to  M  do  begin; 
lor  K:=l  to  I  do 

il  (K<IServar[I])  then  CCI.K]:=K  else  C[I,K] :=ISarver[I] ; 
Pr[old,I,0]:=l; 

Q[I.0]:=0; 
end;  {lor} 

lor  I;=l  to  Paths  do  begin 
lor  K:=l  to  I  do  begin 

il  (K<F JSarver [I])  then  FJCCI.K3;=K  else 
FJCCI.K] :=FJSarverCl] ; 

FJPrtold,I,0]:al; 

FJQCI.O]:=0; 
and;  {lor} 
and;  {lor} 

lor  K:=l  to  I  do  begin 

lor  1:=!  to  3  do  begin 

il  (IServar [I] >=l)  then  RCI,K];=S[I]  else 
il  (IServar [I] =1)  then 
R[I,K] :=SCl]*(l+q[I,K-l])  else  begin 
R[I,K] :=0; 
lor  J: =1  to  K  do 

R[I,K]:=R[I,K]+J*Pr[old,l,J-l]/C[I,J]; 
R[I,K]:=R[I.K]*S[I}; 
and;  {else} 

end;  {lor} 

lor  l:=  1  to  Paths  do  begin 

il  (FJServer[I]>=V[4]'*I)  then  FJRCl.K]  :=FJS[I]  ; 
il  (F JSarver [I]  =  1)  than 

FJRCI.K] :=FJS[I]*(1  +  FJQCl.K-l])  else  begin 
FJR[I,K]:=0; 

For  J : =  1  to  K  do 

FJR[I,K] :=FJR[I.K]  +  J*FJPr [old. I , J-1] /FJC [I , J] ; 
FJR[I,K] :=FJR[I,K]*FJS[I] ; 
and;  {else} 

end;  {lor> 
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for  I : =5  to  H  do  begin 

if  (■Server[I]>=H)  then  RCI.K] ;=S[I]  else 
if  (NServer[I]=l)  then 
RCI.K] :=SCI]*(1+QCI,K-1])  else  begin 
RCI.K] :=0; 
for  J ;=1  to  K  do 

R  Cl .  K] :  =R  Cl ,  K]  +  J*Pr  Cold .  I .  J- 1]  /C  Cl ,  J]  ; 

RCI.K] :=RCI.K]*S Cl]; 
end;  {else} 

end;  {for} 

CalcCyclaTine(Paths .M.V .FJR.FJP .R.CycleTime) ; 

Lambda  CK] ; -K/CycleTime ; 
for  I:=l  to  M  do  begin 

QCI.K]  :=R Cl. K] ♦Lambda CK]*V Cl]; 

UCI.K]  ;=SCl]*LambdaCK]*VCl] ; 
end;  {lor} 

for  I : =  1  to  Paths  do  begin 

F JQ [I , K] : =FJR  Cl . K] sLambda  CK] eFJP  Cl] ‘V  C4] ; 

FJUCl.K] :=FJS Cl] ♦Lambda CK]*FJP Cl] •VC4] ; 
end;  {lor} 

for  I:=l  to  3  do  if  (IServerCl]<l)  then  begin 

for  J:=l  to  K  do  PrCnew.I. J] :=UCl.X]ePr[old,I,J-l]/C[I,J] ; 
PrCnea.I.O]  :=1; 

for  J:=l  to  K  do  Pr Cnes,! ,0] :=PrCnoB,I,0]-PrCnoH.I , J]  ; 
if  (PrCneH,I,0]<0)  then  Error(4.1,K) ; 
end;  {lor} 

PrCold]:*PrCneB]; 

for  I:=l  to  Paths  do  if  (FJServerCl]<VC4]*l)  then  begin 
for  J:=l  to  K  do  FJPrCneo.I.J] :=FJUCl.K]e 
F JPr [old .I,J-1]/FJCCI.J]; 

FJPrCnew.I.O] :=1; 

for  J:=l  to  K  do  FJPrCnew.I.O] :=FJPrCneB,I,0]-FJPr[now.I, J] ; 
if  (FJPr[neB.I.0]<0)  then  Error(6,I,K) ; 
end;  {lor} 

FJPr[old];=FJPrCneB]; 

for  I: =6  to  M  do  if  (IServer [!]<■)  then  begin 

for  J:=l  to  K  do  PrCnes.I , J] :=OCl.K]^PrCold,I,J-l]/C[I. J] ; 
Pr[neH,I.O] :=1; 

for  J:=l  to  K  do  PrCnsB.I.O] ;=PrCnow,I,0]-PrCnew,I, J]  ; 
if  (PrCneH,I,0]<0)  then  Error(4.1.K) ; 
end;  {lor} 

PrCold] ;=PrCneB] ; 
end; 

Briteln( '>  Writing  output 

Report  (I,  M,  Lambda,  irServer,S,V,q,R,U,  Out)  ; 

F JReport ( * . M , Paths , V , Lambda . F JServer . F JS , F J  P , F JQ , F JR . F JU , R , Out ) ; 
close(Dat) ; 
close(Out) ; 

Britain; 

Briteln( ’Program  complete;  output  to  file  KVA. OUT'); 

Britain; 
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and.  {NVA> 


Appendix  C.  Sample  Output 


MVA  PERFORMANCE  REPORT 
Number  of  Customers  (N)  =  24 


St  at 
(i) 

Nuffl  of 

Servrs 

Mean 

Svc  Tm 

Visit 

Ratio 

Thru 

put 

Queue 

Length 

Respns 

Time 

Server 

Util 

1 

99 

0.2500 

1.0000 

3.5432 

0.8858 

0.2500 

0.8858 

2 

99 

2.0000 

0.9500 

3.3661 

6.7322 

2.0000 

6.7322 

3 

2 

1.2500 

0.0500 

0.1772 

0.2239 

1.2637 

0.2215 

4 

1 

1.0000 

0.6675 

2.3651 

0.0000 

0.0000 

2.3651 

S 

8 

1 . 1500 

1.0000 

3.5432 

4.1019 

1.1577 

4.0747 

6 

9 

1.2500 

1.0000 

3.5432 

4.4426 

1.2538 

4.4291 

Average  Queue  Lengths  (including  service) 


N 

i=  1 

i=  2 

is  3 

i*  4 

i*  5 

i»  6 

1 

0.037 

0.285 

0.009 

0.000 

0.172 

0.187 

2 

0.076 

0.569 

0.019 

0.000 

0.345 

0.374 

3 

0.112 

0.854 

0.028 

0,000 

0.517 

0.662 

4 

0.150 

1.138 

0.037 

0.000 

0.689 

0.749 

5 

0.187 

1.423 

0.047 

0.000 

0.861 

0.936 

6 

0.225 

1.708 

0.056 

0.000 

1.034 

1.123 

7 

0.262 

1.992 

0.066 

0.000 

1.206 

1.311 

8 

0.300 

2.276 

0.075 

0.000 

1.378 

1.498 

9 

0.337 

2.661 

0.084 

0.000 

1.550 

1.685 

10 

0.374 

2.845 

0.094 

0.000 

1.722 

1.871 

11 

0.412 

3.128 

0.103 

0.000 

1.894 

2.058 

12 

0.449 

3.412 

0.113 

0.000 

2.066 

2.245 

13 

0.486 

3.695 

0.122 

0.000 

2.236 

2,431 

14 

0.523 

3.977 

0.131 

0.000 

2.407 

2.617 

15 

0.560 

4.259 

0.141 

0.000 

2.578 

2.802 

16 

0.597 

4.540 

0.160 

0.000 

2.749 

2.987 

17 

0.634 

4.820 

0.159 

0.000 

2.919 

3.172 

18 

0.671 

6.099 

0.169 

0.000 

3.088 

3.355 

19 

0.707 

5.376 

0.178 

0.000 

3.258 

3.538 

20 

0.744 

5.662 

0.187 

0.000 

3.427 

3.721 

21 

0.780 

5.926 

0.197 

0.000 

3.596 

3.902 

22 

0.815 

6.197 

0.206 

0.000 

3.764 

4.083 
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23  0.851  6.466  0.215  0.000  3.933  4.263 

24  0.886  6.732  0.224  0.000  4.102  4.443 


Average  Responee  Timee  (including  service) 


N 

i*  1 

i=  2 

i=  3 

i=  4 

i=  5 

1 

0.250 

2.000 

1.250 

0.000 

I.ISO 

2 

0.250 

2.000 

1.250 

0.000 

I.ISO 

3 

0.250 

2.000 

1.250 

0.000 

I.ISO 

4 

0.250 

2.000 

1.250 

0.000 

1.150 

5 

0.250 

2.000 

1.250 

0.000 

1.150 

6 

0.250 

2.000 

1.251 

0.000 

1.150 

7 

0.250 

2.000 

1.251 

0.000 

1.150 

8 

0.250 

2.000 

1.251 

0.000 

1.150 

9 

0.250 

2.000 

1.262 

0.000 

1.150 

10 

0.250 

2.000 

1.252 

0.000 

1.160 

11 

0.250 

2.000 

1.252 

0.000 

1.160 

12 

0.250 

2.000 

1.253 

0.000 

1.150 

13 

0.250 

2.000 

1.254 

0.000 

1.150 

14 

0.250 

2.000 

1.264 

0.000 

1.150 

15 

0,250 

2.000 

1.255 

0.000 

1.150 

16 

0.250 

2.000 

1.266 

0.000 

1.160 

17 

0.250 

2.000 

1.257 

0.000 

1.160 

18 

0.250 

2,000 

1.267 

0.000 

l.lSl 

19 

0.250 

2.000 

1.258 

0.000 

l.lSl 

20 

0.250 

2.000 

1.269 

0.000 

1.152 

21 

0.260 

2.000 

1.260 

0.000 

1.153 

22 

0,250 

2.000 

1.261 

0.000 

1.164 

23 

0.260 

2.000 

1.263 

0.000 

1.166 

24 

0.250 

2.000 

1.264 

0.000 

1.158 
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i=  6 


1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.250 

1.251 

1.251 

1.252 
1.263 
1.254 


4 


Average  Utilizations 


N 

i*  1 

i=  2 

i*  3 

i=  4 

w- 

n 

cn 

i=  6 

1 

0.037 

0.285 

0.009 

0.100 

0.172 

0.187 

2 

0.075 

0.569 

0.019 

0.200 

0.345 

0.374 

3 

0.112 

0.854 

0.028 

0.300 

0.517 

0.562 

4 

0.150 

1.138 

0.037 

0.400 

0.689 

0.749 

5 

0.187 

1.423 

0.047 

0.500 

0.861 

0.936 

6 

0.225 

1.708 

<j  .056 

0.600 

1.034 

1.123 

7 

0.262 

1.992 

0.066 

0.700 

1.206 

1.311 

8 

0.300 

2.276 

0.075 

0.800 

1.378 

1.498 

9 

0.337 

2.561 

0.084 

0.900 

1.550 

1.685 

10 

0.374 

2.845 

0.094 

0.999 

1.722 

1.871 

11 

0.412 

3.128 

0 . 103 

1.099 

1.894 

2.058 

12 

0.449 

3.412 

0.112 

1.199 

2.065 

2.245 

13 

0.486 

3.695 

0.122 

1.298 

2.236 

2.431 

14 

0.523 

3.977 

0.131 

1.397 

2.407 

2.617 

15 

0.560 

4.259 

0.140 

1.496 

2.578 

2.802 

16 

0.597 

4.540 

0.149 

1.595 

2.748 

2.987 

17 

0.634 

4.820 

0.159 

1.693 

2.917 

3.171 

18 

0.671 

5.099 

0.168 

1.791 

3.086 

3.355 

19 

0.707 

5.376 

0.177 

1.889 

3.254 

3.537 

20 

0.744 

5.652 

0.186 

1.986 

3.421 

3.719 

21 

0.780 

5.926 

0.195 

2.082 

3.587 

3.899 

22 

0.815 

6.197 

0.204 

2.177 

3.751 

4.077 

23 

0.851 

6.466 

0.213 

2.272 

3.914 

4.254 

24 

0.886 

6.732 

0.221 

2.365 

4.075 

4.429 
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) 


MVA  PERFORMANCE  REPORT  FOR  FORK-JOIN  QUEUE 
Number  of  Customers  (n)  "24 


St  at 
(i) 

Num  of 

Servrs 

Mean 
Svc  Tm 

Visit 

Ratio 

Thru 

put 

Queue 

Length 

Respns 

Time 

Server 

Util 

1 

3 

2.0000 

0.1742 

0.6173 

1.3159 

2.1317 

1.2346 

2 

4 

2.5000 

0.2043 

0.7237 

1.8838 

2.6029 

1.8093 

3 

4 

3.2500 

0.1328 

0.4707 

1.5639 

3.3229 

1.5296 

4 

3 

3.5000 

0.1215 

0.4305 

1.6862 

3.9173 

1.5066 

S 

5 

4.2500 

0.0988 

0.3500 

1.4927 

4.2644 

1.4877 

6 

5 

4.5000 

0.1435 

0.5085 

2.3384 

4.5986 

2.2882 

Average  Fork- Join  Queue  Lengths  (including  service) 


N 

i=  1 

i=  2 

i=  3 

i=  4 

i=  5 

i-  6 

1 

0.052 

0.076 

0.065 

0.064 

0.063 

0.097 

2 

0.104 

0.153 

0.129 

0.127 

0.126 

0.193 

3 

0.157 

0.229 

0.194 

0.191 

0.189 

0.290 

4 

0.209 

0.306 

0.259 

0.255 

0.252 

0.387 

5 

0.261 

0.382 

0.323 

0.319 

0.314 

0.484 

6 

0.313 

0.459 

0.388 

0.382 

0.377 

0.580 

7 

0.366 

0.535 

0.453 

0.447 

0.440 

0.677 

8 

0.413 

0.612 

0.517 

0.511 

0.503 

0.774 

9 

0.471 

0.689 

0.582 

0.576 

0.566 

0.870 

10 

0.524 

0.765 

0.647 

0.641 

0.629 

0.967 

11 

0.577 

0.842 

0.711 

0.706 

0.691 

1.064 

12 

0.630 

0.919 

0.776 

0.773 

0.754 

1.160 

13 

0.684 

0.996 

0.841 

0.840 

0.817 

1.257 

14 

0.738 

1.074 

0.906 

0.909 

0.879 

1.354 

15 

0.793 

1.152 

0.971 

0.978 

0.941 

1.451 

16 

0.848 

1.230 

1.036 

1.049 

1.004 

1.548 

17 

0.904 

1.309 

1.101 

1.121 

1.066 

1.645 

13 

0.960 

1.388 

1.167 

1.195 

1.128 

1.742 

19 

1.017 

1.468 

1.232 

1.271 

1.189 

1.840 

20 

1.075 

1.550 

1.298 

1.349 

1.251 

1.939 

21 

1.134 

1.631 

1.364 

1.429 

1.312 

2.038 

22 

1.194 

1.714 

1.431 

1.512 

1.373 

2.137 

23 

1.254 

1.798 

1.497 

1.598 

1.433 

2.237 

24 

1.316 

1.884 

1.564 

1.686 

1.493 

2.338 
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Average  Response  Times  (including  service) 


N 

i"  1 

i=  2 

i=  3 

II 

i=  5 

i=  6 

1 

2.000 

2.600 

3.250 

3.500 

4.250 

4.500 

2 

2.000 

2.500 

3.250 

3.500 

4.250 

4.500 

3 

2.000 

2.600 

3.260 

3.500 

4,250 

4.500 

4 

2.000 

2.500 

3.250 

3.500 

4.250 

4.500 

5 

2.000 

2.500 

3.250 

3.501 

4.250 

4.500 

6 

2.001 

2.500 

3.250 

3.503 

4.250 

4.500 

7 

2.002 

2,500 

3.250 

3.606 

4.250 

4.500 

8 

2.003 

2.601 

3.250 

3.510 

4.250 

4.500 

9 

2.005 

2.  jOl 

3.251 

3.515 

4.250 

4.500 

10 

2.007 

2.602 

3.252 

3.523 

4.250 

4.501 

11 

2.010 

2.504 

3.252 

3.532 

4.260 

4.501 

12 

2.014 

2.605 

3.254 

3.543 

1.250 

4.502 

13 

2.018 

2.508 

3.255 

3.557 

4.251 

4.504 

14 

2.023 

2.611 

3.258 

3,573 

4.261 

4.606 

16 

2.029 

2.515 

3.260 

3.591 

4.251 

4.609 

16 

2.036 

2.620 

3.264 

3.613 

4.252 

4.513 

17 

2.044 

2.52t, 

3.268 

3.638 

4.252 

4.518 

18 

2.053 

2,633 

0.273 

3.666 

4.263 

4.524 

19 

2.063 

2.541 

3.279 

3.697 

4.254 

4.532 

20 

2.074 

2.660 

3,286 

3.733 

4.256 

4.541 

21 

2.087 

2.561 

3.293 

3.772 

4.267 

4.552 

22 

2  101 

2.573 

3.302 

3.816 

4.259 

4.566 

23 

2.116 

2.587 

3.312 

CO 

CO 

4.262 

4.681 

24 

2.132 

2.603 

3.323 

3,917 

4.264 

4.699 

(•  (. 


Average  Utilizations 


N 

i=  1 

i=  2 

i=  3 

i=  4 

H- 

II 

cn 

i=  6 

1 

0.052 

0.076 

0.066 

0.064 

0.063 

0.097 

2 

0.104 

0.163 

0.129 

0.127 

0.126 

0.193 

3 

0.157 

0.229 

0.194 

0.191 

0.189 

0.290 

4 

0.209 

0.306 

0.269 

0.255 

0.262 

0.387 

5 

0.261 

0.382 

0.323 

0.313 

0.314 

0.484 

6 

0.313 

0.469 

0.388 

0.382 

0.377 

0.580 

7 

0.366 

0.535 

0.463 

0.446 

0.440 

0.677 

8 

0.417 

0.612 

0.517 

0.609 

0.503 

0.774 

9 

0.470 

0.688 

0.682 

0.573 

0.566 

0.870 

10 

0.622 

0.765 

0.646 

0.637 

0.629 

0.967 

11 

0.674 

0.841 

0.711 

0.700 

0.691 

1.063 

12 

0.626 

0.917 

0.775 

0.764 

0,754 

1.160 

13 

0.678 

0.993 

0.840 

0.827 

0.816 

1.256 

14 

0.729 

1.069 

0.904 

0.890 

0.879 

1.352 

16 

0.781 

1.145 

0.968 

0.963 

0.941 

1.448 

16 

0.833 

1.220 

1.032 

1.016 

1.003 

1.543 

17 

0.884 

1.296 

1.096 

1.079 

1.066 

1.638 

18 

0.936 

1.370 

1.169 

1.141 

1.127 

1.733 

19 

0.986 

1.446 

1.222 

1.203 

1.138 

1.827 

20 

1.037 

1.519 

1.264 

1.265 

1.249 

1.921 

21 

1.087 

1.593 

1.346 

1.326 

1.310 

2.014 

22 

1.137 

1.666 

1.408 

1.387 

1.370 

2.107 

23 

1.186 

1.738 

1.469 

1.447 

1.429 

2.198 

24 

1.235 

1.809 

1.530 

1.607 

1.488 

2.288 
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